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By 
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December  1990 
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Electronically  diabatlc  transitions  arise  from  a 
breakdown  of  the  Born-Oppenheimet  approximation.  Motion  of 
the  nuclei  during  these  breakdowns  can  drive  the  system 
through  changes  in  electronic  configuration. 

Unfortunately,  since  exact  solutions  do  not  exist  for 
systems  with  more  than  a few  degrees  of  freedom, 
approximate  methods  are  needed.  The  elkonal  method 
developed  in  this  thesis  is  a likely  choice  because  of  its 
ability  to  incorporate  any  number  of  potential  energy 
surfaces  in  its  calculations.  The  formalism  has,  In  fact, 
been  developed  and  applied  on  a number  of  occasions  In  the 
diabatlc  representation.  The  present  work  extended  that 
foemallsm  to  the  adiabatic  representation,  and  compared 
application  of  both  sets  of  equations  to  the  sane  model 
systems.  These  included  pairs  of  noncrossing  potentials, 
crossing  potentials  in  repulsive  regions,  and  crossing 
potentials  in  attractive  regions. 

Adiabatic  equations  were  formed  by  applying  a unitary 
transformation  to  the  diabatlc  equations.  The  resulting 
expressions  contained  additional  terms  and  could  not  be 


slspllfied  by  a phase  change  to  the  same  extent  ae  in  the 
diabatic  case.  The  adiabatic  effective  potential  contained 
a nomentun  variable,  P,  that  could  net  be  elialnated. 
Special  measures  were  required  to  form  trajectories  using 
this  momentum-dependent  potential.  P also  appeared  in  the 
denominator  of  the  coefficient  equations,  giving  rise  to 
singularities  in  the  turning  point  region.  Non-Hermitian 
terms  in  the  tine-dependent  equations  made  normalisation  a 
requirement  for  probability  conservation.  Extensive 
calculations  were  performed  of  transition  probabilities 
versus  collision  energies.  Since  both  kinds  of  calculation 
gave  the  sane  results  (as  they  should),  the  less 
complicated  and  computationally  faster  diabatic  equations 
are  preferred.  It  is,  nonetheless,  important  to  have  the 
adiabatic  equations  working  well,  Not  only  can  dynanicists 
learn  from  them,  but  it  is  likely  that  situations  will 
arise  in  which  they  are  useful. 

An  improvement  to  the  usual  electronically  inelastic 
eikonal  calculation  was  also  proposed  and  tested.  It 
imposes  microreversibility  by  combining  the  effective 
potential  for  the  collision  starting  in  the  lower  state  and 
the  effective  potential  for  the  collision  starting  in  the 
upper  state  Into  an  averaged  effective  potential.  The 
improvement  was  dramatic.  Agreement  with  exact 
calculations  was  remarkable.  Because  little  or  no 
additional  computer  time  was  required,  it  was  concluded 
that  eikonal  calculations  are  best  performed  using  an 
averaged  effective  potential  in  the  diabatic  representation. 


CHAFTEIt 


INTRODUCTION 

Eleccconlcally  dlabatic  transitions  arise  fron  a 
breakdown  of  the  Born-Oppenheiner  (Levine,  19E3) 
approsioatlon.  The  basic  situation  is  that  the  electrons 
®ass  is  so  snail  and  their  velocities  so  large  compared  to 
those  of  the  nuclei  that  the  electronic  and  nuclear 
motions  ate  effectively  decoupled.  Variations  in  electron 
velocities  lead  to  situations  in  which  one  or  more 
components  of  the  kinetic  energy  of  the  electrons  goes 
through  tero.  For  a brief  instant  the  electronic  motion 
is  not  so  ouch  faster  than  the  nuclear  motion  Chat  they 
can  be  regarded  as  decoupled.  Any  motion  of  the  nuclei  at 
these  moments  may  drive  Che  system  through  a change  in 
electronic  configuration,  clearly,  the  higher  the  kinetic 
energy  of  the  nuclei,  the  longer  these  moments  will  last, 
and  the  more  chance  there  is  for  a tcansitlon  to  occur. 
Even  low  energy  molecular  collisions,  however,  may  lead  to 
diabatic  transitions  if  the  two  adiabatic  potential  energy 
surfaces  approach  each  other  closely  enough,  some 
chemical  systems  whose  potential  energy  surfaces  intersect 
ace  M » X and  M * Xj  (M  . Li,  Ha,  X,  ...  and  X - F,  Cl, 
Be,  1)  (child,  19791.  In  such  situations,  transformations 
nay  be  applied  that  produce  rigorously  equivalent  sets  of 
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eguations  in  the  dlabstic  and  adiabatic  representations 
(Olson  and  Hicha,  19821. 

Hany  alternatives  to  the  exact  calculation  of 
diabalieally  active  collisions  have  been  devised. 
Prominent  among  these  are  those  that  use  some  sort  of 
senlclassical  approximation.  Reviews  of  seoiclassical 
methods  include  Cross  (1969)  and  Ford  and  wheeler  (1959), 
Reviews  of  diabatlc  interactions  include  Bunker  (1971), 
Child  (1979),  Kicha  (1975),  and  Tully,  (1976). 
Semiclassical  methods  regard  the  motion  of  the  nuclei  as 
predominantly  classical  while  at  the  same  time  maintaining 
the  quantum  superposition  principle  to  include  the  effects 
of  Interferenee  and  tunneling.  Among  the  semielasslcsl 
methods  are  the  MSB  approach  (Messiah,  1961),  the  path 
Integral  formulation  ( Feynmann  and  Hibbs,  1965),  and  the 
complex-valued  trajectory  approach  (Biller  and  George, 
1972).  The  present  work  is  part  of  an  ongoing 
Investigation  into  the  eikonal  method  as  formulated  by 
David  Hicha  in  1983,  improved  upon  in  1986  (Hicha  and 
Swaminathan) , 1967  (Stodden;  Stodden  and  Mlcha),  and  most 
recently  in  1988  (Coheni  Swaminathan,  Stodden,  and  Hicha). 
While  exact  calculations  for  systems  with  more  than  a few 
degrees  of  freedom  ate  not  possible  (Hetman,  1982),  the 
elkonsl  method  is  multidimensional  by  design  and  naturally 
incorporates  any  number  of  surfaces  in  its  calculations. 
Trajectories  are  constructed  that  represent  the  "average" 
state  of  the  system  in  the  Bhrenfest  sense.  These 
trsjectcries  might  well  be  described  as  near-classical  in 


that  they  ace  calculated  icon  Hamilton's  equations  of 
notion  like  the  usual  classical  trajectory,  but  in  a way 
that  relies  on  quantum  coefficients.  The  aethod 


invariably  displays  three  features: 

al  the  wavefunction  used  to  solve  Schcoadlnger' s 
equation  assumes  a forsi  that  is  composed  of  exponential 
and  pre-exponential  factors; 

b)  the  phase  of  that  function  is  associated  with  a 
nechanical  action  integral; 

c)  simplification  is  achieved  with  a short-vavelenqth 
approximation. 


functions  formed  in  this  way  as  pact  of  Che  transition 
integral  for  3-dinensional  elastic  scattering  as  well  as 
in  the  development  of  flux  related  transition  amplitudes 
for  1-dimensional  inelastic  scattering.  The  opening 
discussion  gives  a brief  account  of  the  general  eikonal 
formalism.  Next  is  a description  of  elastic  scattering  as 
recently  treated  by  this  author  (Cohen,  1988).  A point  of 
interest  is  the  use  of  a Jacobian  transformation' to 
facilitate  evaluation  of  the  pre-exponential  factor 
(Xnudson,  Delos,  and  Bloom,  1985;  Delos,  1968;  Stodden, 
1987;  Stodden  and  Hicha,  1967).  The  main  body  follows 
with  a detailed  account  of  Che  eikonal  method  applied  to 
electronically  dlabatic  inelastic  collisions.  Other 


The  major  part  o£  this  voch  was  dedicated  to  applying 
and  conpacing  the  elkonal  approach  in  Che  diabatlc  and 
adiabatic  representations.  (Other  studies  of  this  type 


have  been  done.  See.  for  example. 


1971,  and  sawada  and  Hetiu,  19SS.)  Since  the  eikonal 
method  had  not  previously  been  used  together  with  the 
adiabatic  representation,  the  equations  were  rederlved  and 
new  computational  techniques  developed  for  this  purpose. 

An  improvement  to  the  usual  electronically  inelastic 
eikonal  calculation  was  also  proposed  and  tested. 
Basically,  it  Imposes  the  condition  of  nicroreversibllity 
by  combining  the  effective  potential  for  Che  collision 
starting  in  the  lower  state  and  the  effective  potential 

averaged  effective  potential.  This  implies,  in  essence, 
cunning  the  forward  and  reverse  calculations 


sinultaneously. 


Chapter  II  is  devoted  to  Che  general  formalism  for 
Che  eikonal  method. 

The  first  section  begins  with  the  general  form  of  the 
Schroedlnger  equation  and  ends  with  a more  specialited 
form  involving  nuclear  motion.  The  remainder  of  the  work 
relies  on  this  form. 


The  next  section  starts  with  the  eikonal  form  of  Che 
vavefunction , eliminates  the  exponential  factor  from  the 
nuclear  motion  equation,  and,  by  identifying  Che  phase  as 
the  action  integral,  relates  the  nuclear  motion  equation 
to  Hamilton's  equations  of  notion.  An 


effective 


potential,  related  to  the  quantum  coefficients,  is  formed 
in  the  process.  A rigorous  method  for  specifying  the 
effective  potential  is  given  next.  The  notion  that,  due  to 
conservation  of  flux,  the  coefficients  can  be  related  to 
the  amplitudes  for  transition  probability  is  also 
discussed. 

Chapter  III  is  about  the  eikonal  method  as  applied  in 
the  transition  integral  for  elastic  scattering. 

The  first  section  applies  the  formalism  of  Chapter  II 
to  elastic  scattering. 

The  second  section  starts  with  a general  expression 
for  the  transition  integral  and  develops  an  expression 
that  can,  as  a practical  matter,  be  evaluated.  A Jacobian 

integrated  over  to  be  limited  to  that  spanned  by  a 
representative  bundle  of  classical  trajectories.  The  pre- 
exponential factor  is  also  re-expressed  in  terns  of 
Jacoblans.  The  transition  integral  Is  related  to  the 
differential  cross  section. 

The  third  section  gives  details  on  how  the  Jacobian 
is  evaluated.  It  also  explains  how  the  dimensionality  of 
the  calculation  can  be  reduced  for  an  elastic  collision 
with  a central  potential,  Becauee  no  torque  Is  applied  to 
the  projectile  Its  notion  is  confined  to  a plane. 

The  fourth  section  discusses  computational  aspects 
including  decomposition  of  the  various  vectors, 
utilization  of  the  zero-order  Bessel  function  in 
integrating  over  azimuthal  angle,  «,  and  a linearising 


technique  to  Btabilioe  the  integration 
oscillating  integrand. 


rapidly 


The  fifth  section  displays  a few  sanple  calculations. 

Chapter  iv  develops  the  eikonal  equations  in  the 
diahatic  and  adiabatic  representations. 

The  first  section  shows  how  a unitary  transfornation 
matrix  can  be  applied  to  a diabatie  potential  to  convert 
it  to  an  adiabatic  potential.  Following  the  conversion, 
the  transfornation  angle,  r,  retains  the  inforaation  that 
was  in  the  off-diagonal  elements  of  the  dlabetlo 


The  second  section  applies  and  extends  the  formalism 
of  Chapter  II  to  arrive  at  the  diabetic  equations  for  the 
flux  related  probability  amplitudes.  The  expression  for 
the  diabetic  effective  potential  is  derived.  The 
expression  actually  used  is  simplified  by  the  short- 
wavelength  approximation.  The  specific  form  of  Hamilton's 
equations  are  given.  The  time-dependent  differential 
equations  for  the  coefficients  ace  derived.  Of  special 


Interest  is  a phase  change  that  allows  for  simplification 
of  the  time  dependent  equations. 

The  third  section  applies  the  conversion  presented  in 
the  first  section  to  arrive  at  the  adiebatlc  equations  far 
the  flux  related  probability  amplitudes.  The  expression 
for  the  adiabatic  effective  potential  Is  derived.  Once 
again,  the  expression  actually  used  is  simplified  by  the 
short-wavelength  approximation.  An  unusual  feature  of  the 


the  adiabatic  effective  potential  is  that  it 


contains 


term  in  P,  a raonentum  variable.  Such  aomentum  dependent 
potentials  are  not  unknown  (Goldstein,  1981,  section  1-51 
but  they  do  cooplicate  construction  of  the  classical 
trajectories.  The  specific  fora  of  Hamilton's  equations 
for  the  nuclear  motion  is  qiven  here  also.  Not  only  ace 
they  different  because  of  the  unusual  effective  potential 
but  they  require  inclusion  of  the  spatial  derivatives  of 
the  coefficients.  These  are  entirely  eliminated  in  the 
dlabatic  case.  The  time-dependent  differential  equations 
for  the  coefficients  are  also  derived  In  the  adiabatic 
representation,  but  they  contain  two  different 
singularities.  A very  small  imaginary  increment,  t,  is 
added  to  the  momentum  in  the  denominator  to  avoid  these. 
Finally,  only  part  of  the  phase  change  that  allows  for 
simplification  of  the  time-dependent  diabetic  equations 


can  be  used  with  the  time-dependent  adiabatic  equations. 

The  fourth  section  defines  and  depicts  the  three 
model  potentials  that  are  used  in  the  calculations.  The 
first  one,  called  NC,  has  no  crossing.  The  second  one, 
colled  RC,  has  a ecossing  in  the  repulsive  region.  The 
third  one,  CA,  crosses  in  the  attractive  region. 

chapter  v covers  computational  aspects  of  the 
dlabatic,  adiabatic,  and  exact  calculations.  The  exact 
calculations,  performed  with  the  log-derivative  method, 
are  included  for  comparison  purposes.  Plgures  appear  at 
the  end  of  the  chapter. 


The  second  section  gives  expressions  for  the  first 
derivatives  of  the  transformation  angle,  y,  explains 


how  they  were  cheeked  with  nunetical  derivatives,  and 
shows  how  the  derivatives  of  r can  themselves  be  used  as 
part  of  other  testing  schemes. 

The  third  section  contains  the  espressians  for  the 
adiabatic  potential  elements  and  their  first  derivatives. 
These  derivatives  are  fairly  complicated  and  ace  checked. 

The  fourth  section  explains  how  an  understanding  of 
the  reasonable  behavior  of  classical  trajectories  was  used 
to  piece  together  the  SCE/EP  DIAB  and  ADIAB  ptograns."  It 
goes  on  to  show  that  the  reasonable  behavior  of  the 
effective  potentials  reinforces  the  notion  that  the 
programs  are  working  correctly. 


The  fifth  section  gives  details  on  how  to  use 
LOGDERZV  < the  program  that  does  Che  exact  calculations). 

The  sixth  section  gives  details  on  how  to  use  the 
eikonal  programs. 

Chapter  VI  gives  the  nunetical  results.  Host  of  them 
see  depicted  in  graphs  of  transition  probability  versus 


The  first  section  is  about  the  EXACT  and  SCE/EP  DIAB 
and  ADIAB  (l->2  and  2->l)  calculations.  Features  of 
interest  are  discussed. 


The  second  section  explains  and  demonstrates  Che  use 
of  the  averaged  effective  potential  (AEP)  to  improve  the 
calculations . 


• Self  consistent  Elkonal/Effectlve  Potential  DiABatie  and 
ADIABaCic.  See  APPENDIX  A for  a complete  list  of 
abbreviations . 


The  third  section  gives  the  SCE/AEP  DIAB  and  ADIAB 
results  along  with  some  tables  that  coopare  all  of  the 
calculations.  Some  analysis  is  given.  The  SCE 
calculations  are  compared  to  the  EXACT  ones.  The  tines 
required  for  the  DIAB  and  ADIAB  calculations  ace  also 
conpaced. 

Chapter  VII  is  the  conclusion.  Here  we  discuss  the 
results,  and  point  out  directions  for  future  work. 


CHAPTER  II 

INTRODUCTION  TO  THE  EIKONAL  METHOD 
2-1.  Th»  Schroedlnaer  Equation  tor  wucl««c  Hatlen 


A discussion  of  Che  eikonal  nethcd  may  begin  by 
defining  a generalized  position  variable  Q,  that  appears 
in  the  total  kinetic  energy  operator  K (Hichs,  19B3), 


The  "squaring"  operation  is  actually  a dot  product,  and  o 


States  of  the  whole  system  may  be  described  in  terms  of 
nuclear  positions,  Q,  such  that 


where  IT(Q)>  is  the  total  wavefunction . States  of  the 
electrons  ace  given  within  bracket  natation  so  that 
electron  variables  are  not  explicitly  displayed.  The  ket 
|T(Q)>  represents  <X|T(Q1>  - T(x,  o)  where  X is  the 
collection  of  electronic  variables-  Similarly 


(2-1) 


is  Che  mass  of  each  variable,  after  mass  scaling  them. 


<*IHq|T(Q)>  - Hg(X,  I^J  T(X.  Q| 


IHqI  - b”  ,- 


|T(0)>  - *"(0I 


f Ho  *“(0I  - ^^55  + G^j  ♦"(0}  .(2-51 
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UBlng  Equetlon  (2-5)  twice  in  Equation  (2-21  gives 


V is  the  nucleer  potential  energy  matrix.  A complete 
treatment  might  attempt  to  derive  V from  its  eleetronlc 
structure.  In  what  follows  an  analytical  form  is  usually 
borrowed  from  other  authors'  published  articles. 

Equation  (2-7)  can  be  rewritten  as  De*0,  where  D is  a 
differential  operator.  Kanlpulation  of  this  sinpllfied 
form  allows  one  to  define  a flux  of  probability  that  will 
not  diverge,  and  some  form  of  which  will  satisfy 
conservation.  The  presence  of  in  these  expressions 
leads  to  antl-Hermitian  terms  that  threaten  to  complicate 
a flux  related  description  of  the  dynamics.  A unitary 
transformation  is  introduced, 

I*  (Q)>  - l»“(0)>  L(0)  , (2-8) 

that  is  chosen  to  satisfy  the  imposed  relationship  GeO. 
When  I*  (Q)>  is  the  adiabatic  basis  that  diagonalises  Vg, 
this  choice  for  L leads  to  the  strictly  diabatle 
representation  (Smith,  1969).  The  a’’  are  Che  Born- 
Oppenheimer  wavefunctions  talren  directly  from  electronic 
structure  calculations,  without  G there  is  no  momentum 


coupling  between  different  states.  The  other  diabetic 
representations  that  follow  froB  other  choices  for  L will 
not  be  used  in  this  wort.  The  equation  chat  reoains 
ulcioately  provides  cross-sections  and  transition 
probabilities  for  the  systen  from  the  positions  and 
monenta  of  the  nuclei. 


2-2.  The  flux  Related  Method 


The  wave  function  for  nuclear  motion  can  certainly  be 
separated  into  exponential  and  pre-exponential  factors, 

♦(Ql  - X(OI  exp|lS(OI/hi  . (2-101 

where  X is  a column  multiplied  by  a common  exponential 
factor.  Applying 

(^  Ig)  X exp(lS/K)  - expUS/)tl  |g  * |§  > X U-ID 
twice  in  Equation  (2-9)  gives 


X<0I 


The  derivative  of  the  undeternlned  function,  S(Q),  can  be 
associated  with  a Donentun  function; 

P - 8S/30  (2-lS) 

which  Is  related  to  some  potential,  Uqu(P.O).  as  yet 
undetecoined,  through  the  Hamiltonian  equation 


The  Hamiltonian  itself  is  identioal  to  the  conserved 
(constant)  energy  of  the  system  giving  the  Hanilton-Jacohi 
equation 

55  *Iq’  '’qo'jo’  E (2-n> 

zf  P and  0 are  functions  of  time,  t,  then  dH/dt  - 0 
provided 

8B/3P  . dQ/dt  and  3H/3Q  - -dP/dt  (2-18) 

It  is  clear  from  (2-15)  that  S(Q)  can  be  determined  from  P 
and  g via  the  prescription 


where  the  subscript  refers  to  the  initial  time. 

and  Pj^  are  nuclear  positions  and  momenta  selected 
with  the  collision  partners  separated  by  a great  enough 
distance  to  be  considered  noninteracting.  Specifically, 
the  initial  total  energy,  E,  and  initial  internuclear 
separation,  r^,  ace  pre-selected.  In  the  asymptotic 
region,  changes  in  momentum  for  the  projectile-target 


coordlnaCe,  p-p(r|,  are  considered  ne^lifible,  dp/dC-0  and 
dp/dr-0.  Calculating  gives  p^  through  Equation  (2-16). 
Generally  (though,  as  shall  be  shown,  not  always) 
(dr/dt)j-pj/m. 

To  specify  Uq^^(P,oi  take  Equation  (2-14)  in  the 
sinplified  torn  EX  " 0. 

Left  nultiply  by  X*, 

X'*’PX  ■ 0. 

and  add  this  to  its  own  adjoint, 


Dividing  by  2x'*’X  gives  Equation  (2-17)  with 


j'S(3' 


X w( 


X'30' 


Q)X 


The  specific  forms  of  u for  the  elastic  collision  case,  as 
well  as  for  the  adiabatic  and  diabetic  representations, 
will  be  constructed  below.  This  manner  of  choosing  u 
gives  trajectories  that  represent  the  average  state  of  the 
system  in  the  sense  given  by  the  Ehrenfest  theorem. 

It  remains  to  determine  x,  and  subsequently  x* 
Aotually,  Che  classical  trajectories  are  constructed  from 
time  dependent  equations.  Since  the  coefficients  are 
coupled  to  these  through  the  effective  potential,  a time 
dependent  function,  A(C),  that  is  closely  related  to  Xlr), 
is  used  in  its  place.  The  principle  for  selecting  Che 
initial  values  is  easy  enough,  select  the  initial  state  of 
the  system,  place  all  the  probability  (P)  in  that  state 
For  instance,  if  1 
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lower  state  and  2 the  upper,  and  the  pre-collision 
nolecule  is  in  the  lower,  the  elements  of  A ace  chosen  so 
that  Pj(  lnltial)-1.0  and  Pj  ( initial  1-0 . 0 . Tiioe  dependent 
equations  tor  A ace  basically  transformed  c dependent 
equations  for  X-  As  demonstrated  below  the  diabatic 
equations  are  simpler  because  teems  can  be  eliminated 
thcouqh  changes  in  phase,  The  adiabatic  equations  contain 
singularities  that  prevent  this  kind  of  handling,  so  they 
remain  generally  problematic  and  unwieldy.  Complete 
derivations  will  be  presented  separately  for  the  two 
representations. 


CHAPTER  III 

ELASTIC  COLLISIONS  IN  THREE  DIHENSIQNS 
3-1,  Elestle  Collisions  and  ElEonal  runetlons 


This  chapter  vill  show  hov  the  eikonal  approach  can 
be  applied  to  Che  CransLCion  integral  used  in  conputlng 
cross  sections  Cor  elastic  collisions  (Messiah,  1PS2; 
Bicha,  1963;  Pauly,  1979).  The  discussion  is  sinpllfiod 
by  the  Cact  that  neither  the  target  nor  the  ptojecClle 
have  internal  notions.  They  do,  however,  nova  in  real 
( 3-dlnenslonal } space.  The  kinetic  energy  operator  (see 
Equation  (2-1)),  in  Cartesian  coordinates,  is  represented 


Removing  the  electronic  portion  of  the 
gives  Che  elastic  analog  of  Equation  (2-9), 


(3-2) 

function 


) 


function;  c is 


(3-3 


whece  T(r)  is  the  cadial  wave  function:  c la  the  vee 
cepcesenting  Intecparticle  separation.  Assuming  a 
solution  of  the  fiocn 


T(r)  . x(rl  e‘  , (3_4; 

and  realising  that 

(S  V)  X exp(iS/h)  • expus/h)  ( ^0*78  ) X (3-5) 

(where  the  r has  been  suppressed),  Equation  (3-3)  may  be 
simplified  to 


(k  . 

where  the  "squaring"  of  0«7/i*7S)  refers  to  a dot  product. 
Identifying  the  entire  expression  operating  on  x as  P 
allows  the  total  energy  for  the  elastic  case  to  be 
expressed  as  in  Equation  (2-21).  Pltst,  left  multiply 
Equation  (3-f)  by  x > 


+VX*X  -BX*X  - 0 (3-7) 
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Then,  take  the  adjoint, 

- ?i  I’Vix  - 8 ..  in-,  e , . i|sii  . |i  ,•  f , 

*vx‘x  -BX*X  • 0 . ,3-B) 

Before  combining  the  last  two  equations,  as  in  Equation 
12-20),  it  is  useful  to  consider  S.  The  nature  of  S is, 
as  yet,  undetermined.  Letting  it  satisfy  the  momentum 
equation. 


P - 9S  , ,3_9, 

allows  the  short-wavelength  approximation, 

|d*X/dr^|  <<  |dx/dr  (p/KIl  « |x  (p/)«)^|  , (J-10) 

to  be  applied  to  Equations  (3-7)  and  (3-6).  This 
eliminates  the  terms  containing  V^x  and  VX.  Adding 
Equations  (3-7)  and  (3-8)  eliminates  the  terms  with 
(•9-p)  in  them.  Dividing  the  result  by  2x'x  gives 


Clearly,  p /2m  is  the  kinetic  energy,  so  V(r)  is  the 
effeotive  potential  for  the  elastic  case  after  application 
of  the  short-wavelength  approximation.  Particular  choices 
for  V will  depend  on  the  demands  of  problem  at  hand. 
Anything  from  sn  elaborate  table  meticulously  constructed 
from  experimental  data  to  a simple  Lennard-Jones  guess- 


equation 


Kami Itonianr 


paraaetcically  on  t 
that  trajectories  ma 
equations  of  motion, 


and  equal  to  E.  Both  p and  r depend 
(chosen  here  to  represent  time)  so 
ly  be  constructed  using  Hamilton's 


* ‘l?>r  ■ P''"  ' 

dp/dt  . - (|f)p  - - . (3-13) 


3-2.  EDconal  Functions  and  the  Transition  Integral 


In  many  cases  of  practical  interest,  some  of  which 
ace  Included  in  this  work,  quantum  cross  sections  and 
probabilities  can  be  computed  from  methods  that  rely  on 
conservation  of  flux.  Unfortunately,  when  resonances 
occur,  several  trajectories  with  different  initial 
conditions  (impact  parameters  for  the  elastic  case)  may 
contribute  to  Che  same  result,  while  the  classical  crosi 
section  Is  the  sum  of  such  contributions,  the  quantum 
cross  section  shows  a phase  interference  which  cannot  be 
accounted  for  by  simple  addition. 


The  traneitlon  integral  is  a way  of  calculating 
scattacing  amplitudes  that  specifically  accounts  for  phase 
effects.  It  may  be  represented 


(3-14) 


where  there  is  a plane  wave  associated  with  some  pre- 
selected final  state,  b,  a scattering  wave  associated  with 
the  initial  state  of  the  system,  a,  and  VaV(r)  is  the 
potential  energy.  It  is  the  ability  to  pre-select  the 
initial  and  final  states  that  makes  this  method  so 


appealing.  Dynamics  calculations  can  be  restricted  to 
collisions  with  specified  product  states. 

An  expression  for  the  transition  integral  using 
elkonal  wave  functions  nay  be  derived. 


T . J dr  exp(^  p^.r)  v(c)  xir)  expl^  S(r)J  . (3-15) 

The  first  requirement  for  evaluating  T is  some  way  to 
evaluate  x.  Insetting  Equation  (3-9)  in  Equation  (3-6), 
and  carrying  out  the  operations  indicated  within  the 


brackets 
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The  appearance  of  this  equation  can  be  qreatly  simplified 
with  the  following  relations: 


-|i  . 

® - 7 ■ dc/dt  . d/dr  ■ d/dt  . 
using  these,  Equation  13-16)  becomes 


(3-ni 


iS-20) 


solution 


a simple  first  order  differential  equation  whose 
may  be  expressed  as 


Xir)  • exp(-  i J df  (K^  - iU)j  , (3-21) 

where  tj  refers  to  the  initial  time. 

A high  energy,  short  wavelength,  approximation  allows 
to  be  neglected  when  it  appears  next  to  0,  |K  | <<  |U|. 
This  simplifies  x. 


In  practice,  even  this  simplified  x is  difficult  to 
evaluate.  An  equivalent  expression  that  uses  Jacobians 
has  been  developed  to  avoid  the  problem  (Stodden  and 
Hicha,  1987).  First,  U is  re-expressed. 
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0(t,  . f- 


V is  a three  dinensional  velocity  vector.  The 
transition  integral  of  Equation  {3-11)  is  spatial.  For 
obvious  reasons  the  integration  is  not  carried  out  over 
all  space.  The  useful  information  is  contained  in  the 
region  spanned  by  a representative  bundle  of  classical 
trajectories  generated  by  the  method  of  Section  3-1.  The 
initial  conditions  for  these  trajectories  are  chosen  to 
natch  the  actual  experinental  conditions.  The  integration 
follows  each  classical  trajectory  as  it  proceeds  until  the 
particles  are  out  of  range.  The  time,  t,  the  impact 
parameter,  b,  and  the  asimuthal  angle,  a,  determine  a 
unique  point  along  the  trajectory  that  is  fixed  in  space. 
This  point,  Chen,  can  be  related  to  Che  set  of  cartesian 


coordinates  ( 


y - y{t,b,a) 


where  u is  the  velocity  vector  in  Cartesian  coordinates. 

It  is  worth  noticing  that,  even  though  they  are  determined 
in  two  different  coordinate  systens,  v is  numerically 
equal  to  u^^.  That  is,  the  projection  of  the  velocity 
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veetot  on  the  x-axSs,  a scalar  quantity,  Is  the  sane 
whether  oc  not  that  vector  is  actually  given  in  Cartesian 
coordinates.  In  fact, 

Vy-Uy 

13-25) 

Decooposition  of  V and  v (oc  u)  yields 


where  1,  j,  and  k are  unit  vectors  pointing  along  the  x, 
y,  and  s axis. 

The  three  derivatives  on  the  right  ace  the  terms 
whcse  handling  is  simplified  by  expressing  them  in  terns 
of  Jacobians.  if  used  in  their  present  form,  they  would 
result  in  inpcopec  integrals  buried  inside  other 
integration  loops.  This  would  lead  to  computer  overflow 
and  a host  of  other  mathematical  difficulties.  Taking  the 
first  of  the  three  as  an  example,  we  begin  by  developing 
and  rearranging  the  following  system  of  linear  equations: 

Sx/at  + 3y/St  ^ t Sz/3t  ^ ^ , 

^ * 3y/3b  * 3s/9b  ^ ^ , 


(3-27) 


' • a./..  J . 


3u^/3y  »t.d  atl^/5* 


: 3-311 


/3y  »nd 

Jiilar 

esults. 

1^1^  9y/at  5s/9t 

9x/at  as/at 

7 . V. 

8y/9b  9»/3h 

a^v 

life 

8y/8o  9*/9» 

fefe 

x/Jt  9y/9t  I 
x/9b  8y/9h  | 


»/9.  9y/9< 


But  the  collection  of  nunecatocs  is  ptecisely  3J/tt,  so 


Putting  this  into  U gives 

0(r)  . I (1^]  ln(j)  (3-35) 

A little  Danlpulation  gives 


negative . 
split , 


• exp(-  i I df  (||t]  ln(J)  j - (J/Jj)‘^^*.(3-36) 
1 be  either  positive  or  negative,  the  sguare 
To  locate  Che  coot  on  the  conplex  plane,  J is 


J - |J|  e'  (3-37) 

The  function  arg(J)  is  defined  to  increase  by  n 
every  time  J changes  sign,  x is  now 

)((r)  - (J/Jj)-!^^ 

- expl(-i/21lacg(J)-arg(Jj)n 

(3-38) 

The  exponent  ( acg( J )-arg( ) ! merits  a little 
attention.  Since  J always  equals  at  t^,  the  exponent's 
value  is  zero  at  the  outset  of  every  trajectory.  It  starts 


at  taco  and  increments  by  -in/2  every  time  J changes  sign. 
The  entire  exponential  factor  will  begin  at  1.  go  to  -i. 
-1.  i,  and  then  start  back  at  1,  following  changes  in  the 
sign  of  J.  In  this  way,  the  entire  complex  plane  is 

Getting  back  to  the  transition  integral,  the  various 
expressions  in  the  integrand,  including  the  Jacobian 
transformation,  are  numerically  generated  in  a point-wise 
fashion  along  with  the  trajectories,  with  this 
accomplished,  integration  proceeds  over  time,  impact 
parameter,  and  azimuthal  angle. 


JJJ  dx  dy  ds  f(x,y,z)  - 

The  action  is  rearranged. 


1 - p,.r^  . Jdr. 


* -Pl> 


■ entirely  along  t 
y cosponente.  Th 


sense  because  the  Initial  tine  is  set  before  the  potential 
has  any  influence  on  the  projectile. 

With  X expressed  in  Jacobians  and  the  coordinate 
system  transformed,  the  transition  integral  is  expressed 


I . JJJai  ai.  a.  la^  a|V!  . 

■ -U)  P.,  .J5. 


.(P-Pi)lj  , (3-n> 


. P,.  Pb-  P.  and  pj,  only  ona  c 
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where  Q Is  gensrally  a solid  an9le,  but  in  ths  elastic 
case  is  just  the  scattering  angle. 

3-3.  The  Jacobian  TcansCoraatlon 


A great  sispliftcation,  based  on  a very  sinple 
physical  insight,  occurred  when  the  region  ot  space 
required  by  the  transition  Integral  was  linlted  to  that 
spanned  by  a representative  bundle  o£  trajectories. 

Another  such  insight  also  nakes  things  easier.  A Jacobian 
transformation  from  three-dimensional  Cartesian  to  three- 
dimensional  trajectory  coordinates  has  been  indicated 
(Stodden,  19B7I.  In  the  case  of  an  elastic  collision  with 
a central  potential  no  torque  is  applied  to  the 
projectile.  It  stays  in  whatever  plane  it  starts.  The 
calculation  may  proceed  on  a single  plane  and  the  rest  of 
trajectory  space  recovered  by  formally  rotating  the  plane 
through  o.  Since  every  function  (including  J}  in  the 
transition  integral  will  have  to  be  evaluated  at  every 
point  along  every  trajectory  this  represents  considerable 
savings. 

If  the  plane  of  the  calculation  is  designated  (x’,z) 

X'  - r ain(g)  and  s • r cos{6),  (3-43) 

where  r is  the  internuclear  separation  and  S is  the  angle 
r makes  with  the  z axis.  Since  o is  the  azimuthal  angle 
the  entice  (x'.z)  plane  makes  with  the  (x,z)  plane  when 


placed  in  the  three-dinenslonal  cartesian  coordinate 


A particular  choice  Cor  trajectory  coordinates  uses 
b /2  in  place  oC  b for  the  three-dioensional  case.  This 
causes  the  three-dimensional  sampling  to  be  evenly 
distributed,  since 


(3-45) 


'at  3y/at  aj/3t| 
'3b  9y/3b  3s/3b 
'an  3y/aa  as/a«| 


(3-46) 


Since  zes(a).  az/3a-0.  The  expanded  determinant  becomes 

J - ax/at(-ay/3o  as/3b)l  - ( ay/3t(-3x/3a  3s/3b)l 

* Oz/3t  Ox/3b  ay/3«  - 3x/3a  3y/3b|]  j.  (3-47) 

Since  each  trajectory  is  confined  Co  the  (x'.t)  plane 
am/atBO.  Prom  Equations  (3-35), 


iy/»t  - sin(a). 


9y/3b  ■ 3x’/9b  Binfa), 

»y/»a  - X'  BDS(a).  (3-481 

Using  these  in  Equation  (3-38), 


J - |^|»x'/3t  cos(a)  [-X'  coE(n)  3s/3bl) 

- (8xv»t  sln(ii)  (x'  sin(a)  Ss/Sbl) 

* (as/3t  (5x'/ab  CDS(o)  X'  coe(«) 

♦ X'  sln(«)  Sx'/3b  stn(«)l)j 

Since  sin^(a)  -f  cos^(b)  - 1, 

J - ^ ^3x'/3t  3s/3b  - 3x'/3b  3s/3t]  . 

But  the  expcession  in  Q is  just  > so 

'(;  .v5  a ■ - f -c  a ■ 

Instead  of  evaluating  the  thcee-dlDeneional  Jacobian  in 
ail  three  dinensions,  a two-diaensional  Jacobian  is  used 
and  sinply  multiplied  by  -x'/b. 

Three  of  the  five  factors  in  J'  are  obtained  directly 
from  the  equations  for  motion  in  the  (x',z)  plane.  The 
Kairiiltonian , with  the  momencun  separated  into  components. 


(3-491 


(3-SOI 


(3-51) 


(3-52) 
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The  coupled  equations  of  notion  are  applied  to  each 
coordinate)  f (fax'  or  z),  along  with  its  conjugate 
momentun,  Pj, 

9£/9t  - ( 9H/9pj)  - pj/n 

9Pj/at  - -(9H/9E)  - -OV/at)  (3-531 

Integrating  these  expressions  with  respect  to  tine  gives 
the  factors  z' , 9x'/9t,  and  9s/at  that  appear  in  J'. 

With  the  initial  internuclear  separation  fixed  at 
sone  reasonably  large  value,  the  conjugate  pairs  of 
positions  and  monenta  are  unique  functions  of  the  tine  and 
impact  paraneter  of  a trajectory, 

£ - £(t,bl  pj^  . pj(t,b)  . (3-54) 

The  derivatives  of  Equations  (3-44)  with  respect  to  b are 


Integrating  these  expressions  respect  to  tine  gives  the 
factors  9x'/9b  and  9s/9b  that  appear  in  J*.  J'  and  J car 
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3-4.  CoMDutational  Aspects 

The  problen  of  solving  foe  the  transition  integral  is 
divided  into  two  naicr  tasks:  first,  the  factors  that 
appear  in  the  integrand  are  evaluated;  second,  these 
factors  are  integrated  over  Cine,  impact  parameter,  and 
asinuthal  angle.  The  information  bearing  region  is 
delineated  in  two  steps.  One,  the  internuclear  separation 
at  which  the  trajectories  commence  is  the  distance 
required  for  the  Interaction  potential  to  be  negligible. 
Two,  the  maximum  impact  parameter  has  been  reached  when 
the  flight  of  the  projectile  is  no  longer  Influenced  by 
the  target.  Since  an  elaatic  collision  is  cylindrlcally 
symmetric  the  third  coordinate,  »,  passes  through  3fi0*  in 
all  cases.  A statistically  significant  sample  population 
of  trajectories  is  formed  by  starting  at  sero  Impact 
parameter,  celculating  the  conjugate  poeitions  and 
momenta  for  an  entire  collision.  Caking  small  even 
Increments  In  bV2  and  repeating  the  process  until 
is  reached. 

A packaged  program,  in  this  case  OE,  is  used  to  give 
numerical  solutions  to  the  coupled  differential  equations 
of  motion.  DB  is  a particularly  fine  predictor-corrector 
type  program  written  by  L.  Shampine  and  K.  Gordon 
(Shampine  and  Gordon,  191S).  It  is  well  known  for  its 
reliability  In  a range  of  applications.  All  it  requires 
are  the  initial  function  value  and  derivative  (preferably 
analytical)  for  each  equation  to  be  solved,  a parameter 


i in  Section  3-3. 


. jSt.  #■) . J5..  (i  - ,..J  . 
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The  other  tern  in  the  exponent  of  Equation  (3-32) 
contains  the  only  explicit  reference  to  the  pre-eelected 
final  state.  The  infocDatlon  it  contains  is  critical  to 
the  calculation.  It  also  decostposes  into  tvo  terois. 

C0S(  aH -mVjSin(  9jj)  1 * s nv^  [ l-cos(  9^)  1 . 

(3-57) 

The  magnitude  of  the  momentum  is  conserved  in  an  elastic 
collision,  only  the  scattering  angle,  varies, 


^bx  " sin  ( 0^ ) , 

Pbx  • P P06l6b'  • '3-58) 


By  a convention  that  is  usual  in 
scattering  angles  are  calculated 
Pby"®-  everything  else  is 
a,  integration  over  x reduces  to 


scattering  theory, 
constant  with  respect  to 


I • Jdo  exp{g  lx-  i-mvj)  sin(9^)  cos(«ll)  . (3-59) 


letting 


(l/)() 


(3-«0) 


-J5: 


1(0 . 2,  J5(C|  . 


(3-63) 


T . JJd(bV2)  dt  1(e) 

■ - { i ( J;;-  #■)  • Ip-  - ..-J  - 1 ■ 

• ■ -ili-oxe^ji  j j . 


(3-64; 


The  linearising  technique,  which  nay  be  rendered  as 
nunerically  precise  as  Che  user  wishes,  is  preferable  Co 
the  stationary  phase  approximation.  The  integrand  is 
separated  Into  exponential  and  pre-exponential  cofactors, 

T - Jdt  £(t)  e^  (3-651 


Over  small  increments  of  t, 
be  replaced  by  Che  sero  and 
Taylor  expansions, 


fit)  - f(|  ♦ fjt 


and  6^  are  the  elopes, 


where  i refers  to  inicisl  t 
equal  to  C^  in  Equation  (3-! 


£ and  > are  smooth  enough  to 
first  order  terms  of  their 

(3-66) 


(3-67) 

and  £ to  final  C.  Setting  c 
7)  gives  and 
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Substituting  into  Equation  (3-56), 


- • " ((>.  J;‘  • ‘ • (<i  J;‘ « • ‘ ) ) 

- {( 0 

♦ ( J ^df«jt)  *jt  [ cos(*jt)  + 1 Bin(»jt)l]}. 


in  tho  last  expression,  the  variable  of  integration  w 
changed  fron  t to  4,6  and  the  equality 


e‘“  • eos(u)  * i sin(u)  (3-70) 

was  employed.  The  two  integrals  that  result  can  both  be 
done  analytically. 

The  linearising  technique  integrates  over  £ 
numerically  while  simultaneously  performing  local 
integrations  of  the  exponential  analytically,  integration 
of  a rapidly  oscillating  exponential  is  transformed  into  a 
sequence  of  locally  analytic  integrals  in  the  smoothly 
changing  S. 

After  temporal  integration,  the  integrand  is  complex. 
The  final  integration  over  impact  parameter  is  carried  out 


separately  for  the  real  and  inaginary  parte.  The 
trapezoid  rule  is  used  with  d(b^/2]  approxinated  by 

d(b^/21  . U/2)(b^^j-  ) 

- (l/2l((b.t  ab)^-  b^I  - b^ab  ♦ (i|£i-].(3-71) 

This  ends  the  brief  discussion  of  how  the  theory  ie 
translated  Into  mathematical  expressions  that  can  be 
evaluated  in  the  practical  sense.  Each  of  the  building 
blocks  such  as  the  trajectories,  the  action,  the  Bessel 
function,  or  the  Jacobian  needs  to  be  thoroughly  tested 
both  individually  and  in  various  combinations  before  the 
entire  package  can  be  assembled  in  complete  working  order. 
In  many  ways  producing  a working  program  is  as  difficult 
and  demanding  as  the  theory  itself.  The  interested  reader 
is  referred  to  'Computational  Studies  of  Holecular 
Collisions"  by  this  author  for  more  details  (Cohen,  19681, 

3-5.  Sample  calculations 

The  following  three  figures  depict  sample  eikonal 
calculations  for  elastic  collisions  in  three  dimensions. 

The  first  two  are  for  H.+Hg  at  different  total  energies. 

The  second  one,  at  a higher  energy,  loses  some  of  its 
definition  (the  oscillations  wash  out  somewhat)  but  gives 
a correct  overall  view.  The  third  one  shows  Li.»Hg.  All 
three  were  done  using  a simple  Lenna rd-Jones  potential. 


Figure  3-1.  Quantum  Cross  Section  tor  Hj*Hg  at  R-1,6. 

K is  Che  reduced  energy  (dimensionless). 

The  maximum  error  parameter  (input  to  the  elkonal  program) 

was  10"*. 


An  eikonal  calculation  (solid  line)  compared  to  a partial 
wave  solution  (dashed  line)  provided  by  a program  called 
OPTIK  (1975).  The  results  from  OPTIR  Ire  identical  to 
those  in  Bernstein  (1960), 


DIFFTHENTLfJ.  CROSS  SECTION  (ANG- 
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Figure  3-2.  Quantum  Cto&s  Section  £oc  H2-tK9  at  K*2.6. 

K la  the  reduced  energy  ( dlnenaionless ) . 

The  maxlmuD  error  parameter  (input  to  Che  eikonal  program) 


An  eikonal  calculation  (aolld  line)  compared  to  a partial 
wave  solution  (dashed  line)  provided  by  a program  called 
OPTIX  (197S) . 


(dvaais/ojjv.oHv)  ssojo  i^nosjama 


CEIAPTER  IV 

ELBCTROMICALLY  INELASTIC  SCATTERING  IN  ONE  VARIABLE 
<-l.  Conyertlna  troa  a Dlab»tic  Potantial 


Before  launchin9  into  the  desonstrstion  of  how  the 
elkonal  method  may  be  applied  to  inelastic  scattering  in 
both  the  diabetic  and  adiabatic  representations  it  is 
useful  to  knot*  how  they  intecconvett . The  present  work 
concentrates  on  the  diabetic  to  adiabatic  conversion. 

This  Is  of  particular  interest  because  the  elkonal 
equations  have  already  been  developed  and  applied  in  the 
diabetic  representation.  The  solutions  presented  here  for 
the  adiabatic  case  ace  in  some  degree  extensions  of  these. 
Due  to  the  general  popularity  of  the  diabatie  approach  in 
recent  years  the  literature  contains  a variety  of  diabatie 
potentials  that  can  serve  as  starting  points.  Readers 
interested  in  pursuing  the  reverse  conversion  ace  referred 
to  Olson  and  Hlcha  (Olson  and  Mieha.  1982).  The  unitary 
matrix  Chat  tcanstorms  a diabatie  potential  into  an 
adiabatic  potential  is  described  in  this  section.  When 
the  equation  for  nuclear  motion.  (2-9),  is  written  with  a 
diabatie  potential.  Che  same  matrix  converts  the  equation 
itself  to  the  adiabatic  representation.  This  is  clarified 
in  Section  4-3. 


The  diabetic  potential  in  the  slstple  tvo-state  case 
Is  expressed  as  a two-by-two  matrix  whose  diagonal 
elenents  contain  the  surface  potentials  and  whose  off- 
diagonal  elements  give  the  potential  coupling.  The 
coupling  is  associated  with  the  probability  that  an  atom 
or  molecule  moving  on  one  surface  will  switch  over  to  the 


Another  aspect  of  the  coupling,  that  Is  not  pact  of 
the  potential,  is  momentum  coupling.  Diabatic  potentials 
are  usually  determined  in  such  a way  that  the  momentum 
coupling  can  be  neglected.  Adiabatic  potentials  present 
the  alternate  view,  all  of  the  coupling  is  momentum 
coupling  and  none  of  the  coupling  is  In  the  potential.  By 
definition,  the  off-diagonal  elements  of  the  adiabatic 
potential  are  zero.  Convecslon  from  a diabatic  potential 
to  an  adiabatic  potential  is  effected  by  applying  the 
unitary  transformation  matrix,  B,  containing  the 
transformation  angle,  r,  to  the  diabetic  potential  matrix. 


v®2  . 0 - -vj^sin(r)  eoe(r)  »vJjCob^( 

-vJjSin^lT)  fV^j8in(y)  eos(r) 

- '^2  -''?!>  ==s(r) 


-so- 


Flgutea  4-2,  4-4,  and  4-6  show  the  transformation  angles 


4-2.  gquations  in  the  Dlabatic  Kepresentatlon 

while  more  elaborate  approaohes  are  certainly 
possible,  (surface  scattering,  scattering  in  real  space, 
etc.),  the  present  work  basically  examines  electronically 
inelastic  collisions  in  terns  of  a single  variable,  the 
interparticle  separation  R.  The  development  of  this 
section  generally  follows  Hicha  and  Swaninathan  (1986). 

In  the  diabetic  representation, 

( k Ilf  -6  ) t(RI  ■ 0 . (4-3) 

The  ei)conal  wavefunction  is 

*(R)  - X(R)  expli  S(R)/HJ  . (4-4) 

Ta)cing  the  second  derivative. 


-y  • 


i-5) 
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and  associating  ds/dR  with  the  inoaentusi, 

P - dS/dB  , (4-6) 


- ‘ few  >•  li  ‘ " ■ '“-’I 

An  equivalent  expression  for  this  last  equation  is 
obtained  by  introducing  the  diagonal  matrices  ?r(P,R)  and 
R(B).  H is  introduced  to  facilitate  a phase  change  that 
will  simplify  the  time  dependent  equations.  If  and  K are 
defined 


ll(P,B)  - u(P,R)  + iU(F)  - R(R) 


Equation  (4-?)  becomes 


which  is  still  exact.  Applying  the  short  wavelength 
approximation,  Equation  (3-101,  gives 
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»(F)  • V(B)  + iU(R| 

v(R)  - X^V^X  (4-10) 

The  exact  faco  of  \i“  Is  obtained  as  descEibed  In 
Section  2-2  and  dBoonstrated,  for  the  elastic  case,  in 
Section  3-1.  Identify  the  entire  expression  operating  on 
X in  Equation  (4-7)  as  P.  Left  multiply  that  equation  by 


and  talcs  the  adjoint. 


Adding  these  last  two  equations  and  dividing  by  IX^X  gives 


_ 


-sa- 


Since  R-R(t)  the  function  x<R)  nay  be  expressed  as 
some  function  of  tine.  The  cofactor  A(t)  that  appears  in 
that  function  is  defined  to  satisfy 


XlR(t)!  - exp|ijdf  j Alt)  . (4-14) 

Also  VIR)  • VIRItn  - v^lt),  so  placing  Equation  (4-14)  in 
Equation  (4-10)  gives 


As  pointed  out  in  Section  2-2.  A(t^)-1  in  the  initial 
channel  and  zero  in  all  the  others.  Because  vf(t)  is 
Heroitian. 


d(A*A)/dt  - 0 


A is  conplex. 


a conservation  o 


probability . 


simple  twc-state 


dXj/dt 


dtl/dt  - 0 • <9K/9R)  dR/dt  i-  OH/9F)  dP/dt 
• (9H/3Xjj)  dXj,/dt  * (3H/3Y^)  dY^^/dt 


dR/dt  - 3H/9P  - P/m 
dP/dt  - -3H/aB  • -3\i'^/3B 

dXjj/dt  • 9H/9*n  - 8u'*/aYj| 

dY„/dc  - -8H/ax^  . -au^/aXjj  (4-20) 


In  fact,  the  expressions  for  dx^/dt  and  dY^/dt  in  Equation 
(4-20)  are  identical  to  those  in  Equations  (4-13).  The 
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teader  may  convince  hipself  this  is  so  by  examining  the 
explicit  expcesBlon  foe  v°  given  immediately  following. 
Applying  the  short-wavelength  approximation  and  Equations 


( ''11*1  * ''12*2 

I „d  . . v‘>  A 

'•  '12''2  * “22*2 


• (X^  *fl  1 vjj  * 2 


The  expression  for  -2u^/3R  required  to  generate  the 
trajectories  is  given  by 


-56- 


Notice  that  this  particular  expression  doesn't  contain  the 
derivatives  of  and  3X|^/SR  and  8Y|^/9R.  There  is  a 
tenptation  to  ascribe  this  to  Che  short-vavelenpch 
approximation,  sayinj  terns  that  (for  Instance)  contain 
3X/3R  are  snail  compared  to  terms  in  X,  but  it  really  is 
not  necessary.  Because  Che  expressions  for  dx/dt  and 
dv/dt  used  here.  Equations  (4-lB),  are  identical  to 
Equations  (4-20)  the  inolusion  of  5R„/3R  and  5Y-/3B  in 
iu'^/aR  has  no  net  effect: 


3t^  3r^  * 3t  3r^ 


dXj,  dX^/dt  dY„ 
3r^  " dR/dt  ' 3b^ 


dY„/dt 


aui  ,iuf,  ^ ^ 

3R  'sR  'x,Y  * at  dR/dt 


3t^  dR/dt 


Another  line  of  reasoninq  reinforces  this 
conclusion,  in  the  adiabatic  equations  described  below 
the  expressions  for  dx/dt  and  dY/dt  are  not  equivalent  to 
those  qiven  by  Equations  (4-20).  3X/3R  and  1Y/2R  must  be 

included  in  du/OR  in  that  case  in  order  to  get  results. 
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4-3,  StTuationg  in  the  Adiabatic  Bepcesentation 

The  dlabatic  equation  foe  nuclear  motion,  (d-3),  ie 
conveniently  expressed  as 

( h'^-  e ) . 0 , (4-2dl 


Tcansformation  to  the  adiabatic  representation  is 
accomplished  by  left  multiplying  the  entire  equation  by  b'*' 
and  inserting  the  unit  matrix  in  front  of  if  , 


From  the  definition  of  B,  Equation  (4-lJ,  so 


Equation 


finda 


Inaertlnq  the  first  and  second  derivatives  of  « fron 
Equation  (4-S)  along  with  P-dS/ds,  Equation  {4-6),  gives 
(omitting  from  now  on  the  superscript  a). 


-59- 


;)£}. 

-^•c.  :)s« 

The  first  steps  in  making  this  equation  manageable  are 

a)  apply  the  short-wavelength  approximation  to 
eliminate  d^x/dR^ 

b)  substitute  u for  E - (p^/2») 

cj  group  the  cofactors  of  dx/dR  and  x separately, 


The  term  -i  (h/m)  (dP/dB)  can  he  eliminated  by  a change  in 
phase  similar  to  that  shown  in  Equations  (4-6),  (4-9), 
(4-10),  (4-14),  and  (4-15).  In  principle  one  ought  to  be 
able  to  eliminate  v*  as  well.  In  practice  numerical 
complications  arise  that  prevent  this.  With  the 
expression  multiplying  dx/dR  defined  as  the  c matrix. 


* ‘ X i ( "j  0 ) ai  * 

C can  be  expressed 


Singularities  arise  in  q (and  dx/dB)  near  the  turning 
paint  ICr  as  happens,  |a|  « |0|  or  at  the  turning  point 
when  a - 6 - 0.  The  condition  is  handled  by  adding  a very 
saiall  imaginary  number,  c,  to  P. 


->  q,  - 5 — 5- 

* -(Be(Bj)  +i  10(85)1 

-Be(65)^I«(65)^.''-i2  Rale^)  lo 


-Bo(85>  .10(85)^.12 
-Fe(85)2*lo(85)^«^i2  R*(85)  10(65) 


165)^10(65)^.0^ 


Be(8, 


(see  figures  S-10,  5-U,  i 
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coDplex  let 


The  equation  Coc  dx/dft  becomes 


^ Ib  ) "'"‘Jc'  { ?S  ( 0 -1 


It  is  instructive  to  examine  one  of  the  components 
provided  by  this  equation.  As  an  example,  take 


The  adiabatic  equations  ace  not  simplified  nearly  as 
much  by  a change  in  phase  when  moving  into  the  time  domain 
as  the  diabetic  equations  were.  The  cotactor  of  dx/dR  In 
the  diabatic  equations  was  just  a number  (see  Equation 
(4-9)  ) while  the  cofactor  in  the  adiabatic  case  is  the 
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Bingularities  and  it  becones  impossible  to  find  a 
genuinely  simplifying  phase  change  that  provides 
meaningful  results.  Although  the  adiabatic  equations  that 
remain  are  more  demanding,  the  fact  that  more  of  the 
information  in  the  equation  for  nuclear  motion  is  utilised 

The  adiabatic  effective  potential  u^  is  derived  (as 
in  the  elastic  and  diabatic  cases}  using  Equations  (2-20) 
and  (2-21).  Identify  the  collection  of  terms  operating  on 
X in  Equation  (4-32)  as  f.  Left  multiply  Chat  equation  by 


- «•  ( ", ; ) Si  IS  - ‘ u 

• ^ «•«  - 4 X-  ( “^ ; ) § X 

- ‘ fe  X*«  . |i  x*x  . xVx  - . X*X  - 0 .1. 


i-431 


adjoint 


of  tho  cesultr 
, nl  ,±l,2  . J«i  _*  f “ ‘M  ail  V 

2a  dR  2a  I,  i o J d? 

• i® ' «•  ( ° ‘ ) ii « 

* ^ li  35  " li  ***  " 

i these  last  two  equations  and  divide  by  2x'*'X  to  get 


using  the  short-wavelength  approximation  gives 


<3l'^ 
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Unlike  the  diabetic  effective  potential,  the 
adiabatic  effective  potential  does  not  (after  the  short- 
wavelength  approxioation  is  applied)  reduce  to  a function 
of  R,  but  remains  a function  of  R and  P.  This  has  a 
profound  effect  on  the  equations  of  notion.  In  the  usual 
case  E Is  identified  as  HIR.F)  and,  see  Equation  14-20), 

dR/dt  - 8H/8P  - P/n  (4-47) 

When  u-u(R,P)  however, 

dR/dt  - 8B/SP  - P/n  ♦ 3u/9P  . (4-48) 

The  actual  position  of  the  systen  is  still  given  by 
/dt  (dR/dt)  so  the  klnenatic  momentun  is  m(dR/dt)»P.  What 
is  P?  For  this  particular  set  of  equations  it  suffices  to 
define  P as  ds/dR, 

P a dS/dR  . (4-49) 

it  is  now  convenient  to  give  u®,  dR/dt,  and  dp/dt 
explicitly: 


:;ij  { [ 


. 2 H 


. i2 


Despite  the  additional  tern  in  dR/dt,  it  is  still 
fairly  compact, 

da/dt  - 8H/SP  - p/n  + 8V/9P 

- (P/m)  + 2 (tin2-  Jjtii)  . (4-52) 

The  adiabatic  expression  for  dP/dt  contains  terns  with 
df|^/dR  and  dh^/dR  that  do  not  cancel  each  other  (see 
Equations  (4-18)  and  (4-20)  and  the  surrounding  text)  and 
it  is  fairly  lengthy, 
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dP/dt  - -JB/SB  - -S«®/9R 


Since  in^  it  is  tempting  to  neglect  terms  in 

dl./dB  and  dh^/dB  using  the  short-wavelength 
approximation.  Actual  calculations  reveal  Chat  this  is 
inaccurate.  Apparently  the  short  wavelength  approximation 
doesn't  apply  In  this  situation. 

Changing  to  the  t variable,  since  the  adiabatic 
equations  weren't  as  simplified  by  changing  the  phase,  the 
adiabatic  analogue  of  Equation  (4-15)  is 
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Since  H is  not  HecmiCien  a"*’a  is  not  auconatically 
conserved , 


Conservation  of  probability  is  achieved  through 
nornalizatlon . The  normalised  components,  X'  and 


so  the  total  normalized  probability  is  one. 


(«-57) 


(4-56) 


It  is  clear  from  the  material  presented  in  this 
section  that  there  are  a number  of  drawbacks  to  applying 
the  adiabatic  equations.  First,  there  ace  the 
singularities.  Attempting  to  deal  with  the  singularities 
leads  to  other  complications.  Not  only  is  there  a mote 
complicated  set  of  equations  to  actually  solve,  but  there 
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is  c to  deal  with.  It  must  somehow  be  demonstrated  that 
the  selected  value  of  e is  reasonable.  No  matter  how  this 
is  done  it  requires  additional  work,  coefficients  that 
require  normalisation  to  conserve  probability  further 
complicate  natters.  The  complicated  adiabatic  equations 
are  clearly  less  hardy  and  resilient  than  the  simple 
symmetric  diabatic  ones  and  should  almost  always  be 
avoided.  On  the  other  hand  being  able  to  solve  the 
adiabatic  equations  represents  another  capability/  another 
tool  as  it  were,  in  the  dynamicists'  methodology. 

Somewhere  in  the  great  variety  of  situations,  the 
adiabatic  representation  is  bound  to  find  appropriate 
application. 

4-4.  The  Model  Potentials 


Three  model  matrix  potentials  were  selected  to 
represent  the  general  situation:  a noncrossing  potential 
(NO,  a potential  whose  diagonal  elements  cross  in  the 
repulsive  region  (CR),  and  a potential  whose  diagonal 
elements  cross  in  the  attractive  region  (CA).  The 
noncrossing  potential  was  formed  by  adding  the  terns  in 
the  upper  curve  of  CR  rather  than  taking  their  difference. 
The  resulting  curves  ace  separated  enough  to  give 
transition  probabilities  that  could  be  too  small  to  serve 
as  a useful  example.  To  compensate,  the  coupling 
parameter  a^2  chosen  equal  to  0.05.  The  NC  potential 
still  remains  something  of  a challenge  to  Che  sensitivity 


of  Che  calculations.  An  exaople  of  a nonccosslng 
potential  Is  He*  * H IZygelman  et  al.,  1989).  Presumably 
there  are  other  A * B systems  vith  this  same  kind  of 


potential.  The  HC  case  has  been  discussed  in  the 
literature  under  Che  name  Demkov's  model. 

The  CS  potential  itself  was  taken  from  Currier  and 
Herman  (Currier  and  Herman,  1985).  Ronald  McCarroll,  in 
paper  on  ion  atom  collisions  (Kccarroll,  1980),  depicts 
several  potentials  that  cross  in  Che  repulsive  region. 


The  N + H and  CH*  potential  energy  curves  both  have 
such  crossings  (see  Figs.  7-2  and  9-9  in  McCatroll's 

repulsive  region. 

The  CA  potential  is  from  Savada  and  Hetiu  (Sawada  and 
Hetiu,  1986).  This  is  the  type  of  potential  (for  simple 
systems  anyway)  most  frequently  repotted  on  in  the 
literature.  The  entice  class  of  H X (where  H is  a metal 
and  X is  a halogen)  interactions  has  a characteristic 
potential  Chat  crosses  in  the  attractive  region.  Similar 
crossings  arise  in  atom-molecule  systems  such  as  K X, 
and  n + no,  Oj,  SFj,  . . . and  many  others  (Baede,  1975). 

Even  relatively  simple  scattering  events  ocasionally 
require  a multitude  of  surfaces  Co  describe  them 
adequateiy.  The  example  shown  next  (Olson  and  Garrison, 
1985)  has  fewer  porCential  curves  than  one  would  like  to 
handle  comfortably. 


Adiabatic  potentials  foe  tSe  sodiuB  atoo  intecaeting 
the  Mllioi  surface.  'Olson  and  oaccison,  198S;. 


The  following  tables  give  all  the  relevant 
information  for  each  model  potential,  xn  addition  to 
tables  ace  figures  giving  diabatic  {solid  lines)  and 
adiabatic  idashed  lines)  pictures  of  each  of  them. 
Accompanying  the  potential  figures  are  graphs  of  the 


transformation  angles 
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TABLE  4-2.  POTENTIAL  WITH  CROSSING 
IN  REPULSIVE  REGION  (CR) 
equations: 

V^^-  A^^  e»p(-B  R) 

V - A exp(-o  R)  - a exp(-6  R)  ♦ V 
22  22  22  22  22  20 

Vi2  ■ ^ ) 

paraaeters: 


a ■ 3.0  Bohc  ; B - 1.0  H 
22  22 

12  '**  ' '’20  ’ 


TABLE  4-3.  POTENTIAL  WITH  CROSSING 
IN  ATTRACTIVE  REGION  (CA) 


( exp(-2o(B-B  ) t 2 exp[-a(H-R  ))  | - (e/2| 


- 0 I exp[-2«(B-B  ) - 2 expl-oiH-H  H ) ♦ (e/2) 


paranetecB : 


•02S  K.  i c • .OOS  H.  ! B - 0.4 


• atomic  units  (H.«  Hactceel 
••  unitless 


SUBTpBj//. 
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CHAPTER 


COMPUTATIONAL  ASPECTS — 
ADDITIONAL  DETAILS  OF  THE  CALCULATION 

5-1.  Introduction 


This  chapter  is  about  software  development  end 
implementation.  The  second  and  third  sections  give 
analytical  expressions  for  the  various  derivatives 
required  in  the  SCE  programs.  There  are  the  first  two 
derivatives  of  y,  required  by  the  adiabatic  programs.  The 
first  two  derivatives  of  the  diabetic  potential  elements 
are  needed  to  get  the  derivatives  of  y.  The  first 
derivatives  of  the  diabatic  and  adiabatic  potential 
elements  are  needed  as  well  for  SCE  DIAB  and  SCE  ADIAB 
respectively.  An  alternative  to  taking  the  analytical 
derivative  directly  is  to  spline  fit  the  function  and  use 
the  coefficients  of  the  spline  fit  function  to  form 
numerical  derivatives.  This  is  a natural  way  to  check  the 
analytical  derivatives.  Another  way  is  to  use  the 
derivatives  as  part  of  some  function  that  can  be 
integrated  numerically  by  the  differential  equation 
solving  program  (in  this  case  DE).  This  has  the  advantage 
of  simultaneously  checking  the  solver  itself.  Both 
methods  were  used. 
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The  foutth  section  first  explains  how  insights  about 
the  behavior  of  classical  trajectories  (energy 
coneervatioRr  continuity,  etc.)  could  be  used  to  develop 
the  adiabatic  oikonal  prograos  in  a piecemeal  fashion. 
Starting  with  constant  coefficients  and  a simple  malieBhift 
effective  potential,  more  complex  and  realistic 
expressions  were  added  bit  by  bit  until  the  formulae 
developed  In  Section  4-3  were  entirely  included  in  the 
program  and  producing  reasonable  loolcing  trajectories. 

The  fourth  section  goes  on  to  show  that  the  various 
effective  potentials  behaved  reasonably,  thus  reinforcing 
the  validity  of  the  program.  Since  the  diabetic 
equations,  developed  at  an  earlier  date,  were  already 
known  to  produce  reasonable  results,  much  consideration 
was  given  throughout  to  similarities  in  behavior 
demonstrated  by  the  adiabatic  equations. 

The  fifth  and  sixth  sections  ace  devoted  to 
operational  details  of  the  exact  (LOGDBRIV)  and  eikonal 
(SCE)  programs  respectively.  They  ace  pretty  direct  and 
self-explanatory. 


Derivatives  of  y 


In  addition  to  the  transformation  angle  r,  Equation 
(4-2),  the  adiabatic  equations  for  the  coefficients 
require  the  first  and  second  derivatives  of  y with  respect 
to  the  intsrnucleac  separation  R (see  Equation  (4-55)  ). 


since  T is  a function  of  the  diabatic  potential  elenents 
''ll'  ''22'  '^2'  fleet  and  second  derivatives  of 

these  elenents  go  into  the  first  and  second  derivatives  of 
f.  A sinple  hierarchy  of  derivatives  can  be  constructed 
for  the  MC  potential: 


v^l 


Tj2  ■ tanhl  (R  - Rjj ) j 


T (1 


(5-1) 
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A convenient  expression  for  the  first  derivative  of  t was 
found  in  an  article  by  Evans,  Cohen,  and  Lane  (1971). 
Taking  the  derivative  of  this  gave  d^r/dB^  . 


(5-2) 
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Th«se  expressions  are  certainly  conpllcated  enough  to 
merit  some  cheeking  and  testing.  The  first  step  in  the 
testing  Is  accomplished  by  evaluating  the  function  along 
with  its  analytical  derivatives  over  the  Interesting  range 
of  a.  The  function  values,  along  with  the  values  of  R, 
ace  input  to  an  international  Hathematics  and  Statistics 
Libcary  subroutine.  ICSBVU  (ICSCCO,  1978),  that  treats 
then  as  data  points  and  forms  a cubic-spline  fit  to  the 
original  function.  The  coefficients  of  the  spline  fit 
function  are  used  to  provide  numerical  derivatives  to  a 
high  degree  of  accuracy  ( IHSL  subroutine  QCSEVU  (DCSEVU. 
1961}  ],  The  analytical  derivatives  are  compared  to 
these.  The  testing  follows  the  hierarchy  presented  above: 
first  each  potential  element  and  its  derivatives  are 
tested,  then  v and  its  derivatives.  Figures  5-1,  5-2,  and 
5-3  show  the  final  results  of  this  testing  for  the  NC 
potential.  The  derivatives  were  tested  for  all  three 
potentials  until  analytical  and  numerical  values  were  as 
graphically  indistinguishable  as  in  Che  NC  figures.  The 
equations  for  the  CB  potential  ace  the  same  as  tor  the  NC 
except  that  is  used  in  place  of  -Bjj  and  the  coupling 
parameter  is  0.01  instead  of  0.05. 

Following  are  the  derivatives  for  the  CA  potential: 

?i  a Uj  exp[-2«(B  - B(|ll  : Tj  ■ Op  2 exp[-o  (B  - Hp)l 

* Information  on  insL  software  nay  be  obtained  at  7500 
Bellaire  Blvd.,  Houston  TX  77036. 


''tl 

dvjj  dTj  dTj 

35  " 35“  * 35” 


The  analytical  derivatives  of  r vete  also  used  in  an 
entirely  different  kind  of  testing.  In  addition  to 
evaluating  y directly,  Equation  (4-2),  there  is 
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If  everything  Is  working  as  It  should  the  differential 
equation  solver  (DE)  will  integrate  dr/dB  dR/dt  over  tine 
and  obtain  the  sane  value  of  t as  gotten  from  Equation 
(4-11.  If  the  expressions  for  either  dr/dR  or  dB/dt  are 
sonehow  incorrect  or  OB  is  not  operating  correctly,  the 
T's  will  not  match.  This  is  a very  sophisticated  test 
that  can  help  locate  a number  of  different  problems.  It 
is  even  more  powerful  when  coupled  to  a tost  for  dr/dR. 
dr/dR  is  evaluated  analytically  and  also  from 

If,  for  instance,  the  test  for  r succeeded  while  the  test 
for  dr/dR  failed,  it  would  be  evident  that  d^r/dR*  was 
incorporated  incorrectly  into  the  program,  figure  5-4 
shows  typical  results  from  this  type  of  testing.  Usually 
graphs  would  not  be  necessary,  one  would  simply  monitor 
the  numbers  along  with  other  intermediate  results  as  the 
calculation  proceeded.  This  testing  does  require  the 
solving  of  two  additional  differential  equations,  once 
the  program  tests  correctly  they  ace  removed  to  save 
computer  time. 


DerivativB 


the  Adiabatic  Potential  Elaagnta 


Fron  Equation  (4-1)  come  the  adiabatic  potential 
elements , 

Vjj  ■ cos(y)  I eo6(rl  - sin(T)  ) 

-Bln(T)  I -sln(Tl  Wjj  ♦ eosCyl  ) ; 

* sin(T)  I slnlvl  vjj  -f  eos(r)  Vjj  1 

♦eoslYi  t sin(v)  + coa(r)  v|j  I (S- 

The  time  derivative  of  P,  Equation  (4-22),  reguices  the 
first  derivatives  of  these  eiements  with  respect  to  s, 

jgii  - -sin(Y)  [ cos(y)  Vjj  - sin(r)  vjj  ) 

*eos(Y)  I -sln(Y)  vfi  * oos(Y)  ^ 

- C08(y)  ^ vjj  - sin(Y)  1 

-cos(y)  I -sln(Y)  v|j  Y cosIy)  vfj  1 

jj^d 

-sin(Y)  i -cosIy)  Ib  Vjj  - Sin(Y) 

- sin( y)  ^ v^2  * 
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as^  - co3(tI  3^  I ain(T)  + cosCtI  1 

. . dy? 

♦ sin(y)  1 C06(t)  t sitilr)  gj — 

-Bin(r)  gj  I cos(y|  ♦ sintr)  V^j  1 


+cos(r)  ! -sin(T)  Vjj  + cob(t)  gg^ 

* coa(r)  g^  vjj  * sin(yl  1 (5-7) 

Hhile  thaae  expressions  are  fairly  easy  to  derive,  they 
are  also  lengthy  and  were  tested  by  comparison  with 
numerical  derivatives  aB  in  the  first  part  of  the  previous 
section,  rigure  5-5  Bhows  the  resulte  of  theee  teets. 

5-i.  Reasonable  Behavtor- 
Classical  TraiectorieB  and  Ettectlve  Potentials 

In  addition  to  exact  tests,  a large  part  of  being 
able  to  develop  and  work  with  programs  like  these  cones 
from  having  some  idea  of  what  constitutes  reasonable 
behavior.  The  differential  equations  for  the  adiabatic 
coefficients  are  extremely  complicated  and  closely 
intertwined  so  that  any  small  error  in  any  of  then  could 
easily  destroy  Che  whole  calculation.  Trying  to  get  these 


right  whils  simultaneously  developing  classical 
trajectories  from  momentum-dependent  potentials  for  the 
first  time  {for  this  author!  proved  to  be  impossible.  The 
problem  was  broken  into  small  pieces.  As  each  piece  was 
handled  in  what  appeared  to  be  a reasonable  way  the  next 
piece  was  added. 

The  first  step  was  to  calculate  trajectories  from  a 
simple  momentum-dependent  effective  potential  with 
constant  coefficients.  As  no  relation  to  any  physical 
reality  was  intended  any  improvised  potential  would  serve. 
The  actual  choice  was 


When  this  calculation  developed  to  the  point  where  the 
total  energy  was  conserved  and  the  trajectory  was 
continuous,  s didn't  increase  or  decrease  haphasardly,  the 
next  step  was  taken.  A single  simple  non-constant 
coefficient  was  added.  Again,  there  was  no  need  to  be 
realistic,  any  simple  expression  for  the  coefficient  could 
be  used,  when  this  new  trajectory  demonstrated  reasonable 

added.  Things  were  gradually  built  up  in  this  manner 


until,  after  a week  or  so,  all  the  actual  (derived) 
expressions  were  in  place  producing  what  appeared  to  be 
reasonable  classical  trajectories,  while  there  is  no 
reason  to  reproduce  this  process  here,  understanding  the 
approach  may  be  useful  to  the  reader  who  intends  to 
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pertoco  these  calculations.  Fi9uree  S-6  and  S-7  conpace 
the  positions  and  momenta  for  a typical  trajectory 
obtained  with  SCE/AB  OIAB  and  SCB/AE  ADIAB.  The  agreement 
lends  credibility  to  the  adiabatic  calculations  since 
dlabatic  calculations  are  already  fairly  well  understood. 

Another  clue  that  the  calculations  are  behaving 
reasonably  can  be  gotten  by  examining  the  effective 
potentials.  The  non-averaged  effective  potential 
generally  (at  Least  in  the  asymptotic  regions)  reflects 
the  transition  properties  of  the  system.  Suppose,  for 
example,  one  were  examining  a l->2  transition.  As  long  as 
the  transition  probability  remained  zero  the  effective 
potential  would  remain  equal  to  the  ground  state  potential 
If  the  transition  probability  were  to  become  (let  us 
say)  0.2,  the  effective  potential  would  relocate 

about  two-tenths  of  the  way  from  to  If  one  were 

looking  at  a 2->l  transition  it  would  locate  two-tenths  of 
the  way  down  from  the  upper  surface  rather  than  two-tenths 
of  the  way  up  from  the  lover.  Reasonable  behavior  for  the 
averaged  effective  potential  is  quite  different.  Due  to 
its  symmetric  nature  one  can  expect  it  to  remain  about 
half  way  between  the  upper  and  lower  surfaces  at  all 
times.  Its  location  does  not  reflect  the  transition 
probabilities.  These  can  only  be  gotten  from  Che 
amplitudes . 

Figure  5-8  shows  the  three  effective  potentials  for 
three  SCE  DIAB  calculations  at  a total  energy  of  1.75 
Hartree.  The  bottom  one  is  for  an  SCE/EP  l->2 


calculation.  Consulting  Plguce  6-2  ceveals  that  (fron 
the  l->2  calculation)  at  this  energy  is  extremely  small. 
Accordingly,  Che  eftective  potential  stays  close  to 
The  top  curve  is  for  an  SCE/BP  2->l  calculation. 

Consulting  figure  6-2  again  reveals  Chat  P^2  the 

2->l  calculation)  at  the  saoe  energy  is  close  to  0.2, 
P22SO.2.  The  effective  potential  in  this  case  goes  in  at 
Che  V22  level  but  comes  out  about  two-tenthe  of  the  way 
down.  Finally,  Che  SCE/ASP  DIAB  averaged  effective 
potential  goes  in  and  comes  out  precisely  half  way  between 
Che  upper  and  lower  surfaces.  Figure  5-9  shows  the  same 
information  Cor  the  5CE/AEF  ADIAB  calculation.  (The 
appropriate  transition  probabilities  nay  be  obtained  Cron 
Figure  6-3.)  The  remarkable  agreement  with  Che  diebacic 
calculations  once  again  reinforces  Che  credibility  of  the 
adiabatic  calculations. 

On  the  whole,  the  effective  potentials  do  behave 
reasonably  for  both  cases.  The  only  real  difference  is 
Chet  the  outbound  adiabatic  averaged  effective  potential 
is  not  generally  precisely  half  way  between  the  two 
surfaces.  This  is  entirely  an  effect  of  passing  through 
the  singularity  at  the  turning  point.  To  prove  this, 
increase  e to  the  point  where  the  singularity  isn't  sensed 
by  the  equetions,  and  the  averaged  effective  potential 
will  remain  at  the  mean  value  (even  though  the 
coefficients  themselves  show  the  normal  range  of 
oscillation).  Fortunately,  because  of  the  way  the 


coef flcienCs 


selected. 


acceptable  valuea  of  the 
erstwhile  Inperfeetlon  has  no  effect  on  the  calculation. 

5-5.  Exact  Calculations — 


This  work  usee  a nuDber  of  results  foe  coapacison 
purposes  that  are  produced  by  calculations  from  exact 
methods.  These  methods  do  not  in  general  use  the  tine 
dependent  equations.  The  particular  one  selected  for  this 
project,  the  log  derivative  method,  is  well  understood  and 
widely  used  and  will  not  be  explained  here,  interested 
readers  are  referred  to  Lawrence  Relyea's  doctoral 
dissertation  (Relyea,  1980),  an  article  by  Johnson 
(Johnson.  1973),  along  with  the  host  of  other  references 
found  in  current  literature.  The  progran  itself, 

LOGDEBIV,  was  written  by  Pr.  Belyea  in  1977.  It  requires 
only  a few  input  parameters  and  is  extremely  easy  to  use. 
The  following  tables  summarize  input  for  the  three  model 
potentials  (an  explanation  follows  Table  5-1).  In  fact,  a 
wide  variety  of  forms  for  the  output  also  require  some 
selection  at  the  input  level. 
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Table  5-1.  Input  for  LOGDEBIV  - NC  and  RC  Potentiala 


VI  - V2  . V3  - VS  - V6  • V7  . V8  - V9  - no  - D.O 
V4  . 2.0  Bohr 

EPS(l)  - 0.0  Hacttee  ; EPS(2)  • 0.1  Hartree 
ELOSTdl  . 0.0  Hartree 
H(l)  - H(2I  • 0.001  Bohr 


NH  . 1 i NBLOST  • 1 i ISPY  - I ; IVTYPE  - 17 

B is  the  total  energy.  It  is  left  blank  here  because 
different  ranges  were  used  for  different  calculations. 

RBDHAS  is  the  reduced  mass  of  the  systeD. 

YO(I)  is  the  initial  value  of  the  wave  function. 

LHIN  and  LHAX  refer  to  the  mininun  and  naxiaum  values 
of  the  angular  raoaentum  guantuo  number,  presumed  for  these 
calculations  to  be  sero. 

VI  to  VIO  are  potential  parameters.  Since  the 
present  potentials  were  added  to  the  existing  program 
these  weren't  needed.  V4  is,  however,  a characteristic 
length  used  by  the  program  Itself  for  locating  turning 
points.  Generally  the  location  of  a potential  minimum  is 
used  for  V4.  This  needn't  be  exact. 


SPSCl]  gives  the  asyuptotic  energy  of  the  Ith 
potential  surface.  In  the  usual  case  the  lover  surface 
at  the  teto  of  energy,  EPS(I)  - 0.0  hartcee.  In  that 
EFS(2)  gives  the  separation  of  the  two  surfaces. 

GLOST(I)  represents  Che  energy  lost  with  the 
electron.  No  electrons  ace  lost  in  these  calculations. 

H(I)  is  the  stepsize  used  for  the  calculation  on  t 
Ith  surface,  The  value  given,  0.001  bohc,  snail  as  it 
seem,  is  the  largest  value  at  which  the  calculations 
consistently  converged.  This  is  important  because  it 
increases  the  computer  time  required  by  this  method. 

xniN  and  XKAX  represent  the  minimum  and  oaximun 
values  of  the  internuclear  separation  used  in  the 
calculation.  Probably  a little  smaller  value  of  xnAX 
could  have  been  used,  but  not  much. 

NCHN  is  the  number  of  open  channels. 

NH  is  not  set  to  one  only  if  several  values  of  the 
stepsize  ace  used;  NELOST  is  not  set  to  one  only  if 
electrons  ace  lost;  ISFY  is  not  set  to  one  only  if  the 
matrix  is  calculated  from  Che  R-oatrix  rather  chan 
directly  from  the  exact  solution.  Everything  else  stai 
at  zero  except  IVTYPE,  which  didn't  apply  because  new 
were  inserted. 


po tentials 


5-6.  Eikonal  Calculations — 


i eikonal  calculations  covered  in  the  present  work 
the  most  pact  exercises  in  solving  sets  of  coupled 


differential  equations, 
prepackaged  progran  DE  1 
advantage  in  using  DB  is 

computer  time.  Any  of  t 


As  previously  stated,  the 
used  for  this  purpose.  The 
its  ability  to  vary  Che  siae  of 
saves  a considerable  amount  of 
e eikonal  programs  themselves 
>2  or  2->l,  SCE/AEP  DIAB  or 


ADlAfi,  NC,  CR,  or  CA)  can  be  divided,  as  far  as  the  user 
is  concerned,  into  three  major  pacts.  There  is  PROGRAM 
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MAIN,  which  contains  the  initial  conditions  for  all  the 
difCerential  equations  along  with  a few  operating 
paranetecs  required  by  DE,  the  calling  statenent  to 
Subroutine  DE,  and  provisions  for  putting  the  output  into 
a desired  torn  (graph,  table,  etc.  ).  There  is  Subroutine 
F that  contains  the  temporal  derivatives,  analytical  in 
the  present  case,  for  each  equation  to  be  solved  by  DE. 
There  is  Subroutine  DE  itself  that,  with  the  single 
exception  discussed  below  (Rl),  needn't  concern  the 
present  user,  in  addition  to  the  inputs  to  DE  one  other 
input,  the  value  of  e,  is  required  for  the  adiabatic 
calculations. 

Those  interested  in  looking  at  an  eikonal  program  are 
referred  to  APPENDIX  s.  which  contains  a fully  commented 
printout  of  the  original  portions  (not  DE)  of  the  code  for 
SCE/AEP  DIAB.  The  discussion  that  follows  is  more 

Tables  S-3  and  5-4  show  physically  related  input, 
some  explanatory  material  follows  Table  S-3.  A discussion 
of  program  related  input  follows  Table  5-5.  Finally, 
there  is  a discussion  of  e. 

Table  5-3.  Input  to  DE  - HC  and  RC  Potentials 


RH  - 1636.15  au 


All  of  the  pcesenC  SCE  calculations  ace  written  and 
pecfomed  entirely  in  atomic  units.  E is  the  total  energy 
in  Hartcees,  left  blank  here  (as  In  the  last  seccioni 
because  its  range  varies  with  the  calculation. 

RI  is  the  initial  value  of  the  internuclear 
separation,  chosen  to  be  in  the  asynptotlc  region  of  both 
potential  surfaces  and  well  past  any  crossing  (or  avoided 
crossing) . Rl  also  furnishes  the  cutoff  value  for  the 
calculation.  This  Is  unusual  in  the  way  bE  is  handled. 

In  the  customary  use  of  DE  a maxinuo  cutoff  time  is  passed 
through  the  calling  statement  along  with  the  rest  of  the 
input.  This  value  is  fixed  and  cannot  respond  to  the 
calculation.  In  the  present  case  this  parameter  (called 


e large  dummy  value  t 


reached.  A flag,  NEND,  is  set  to  sero.  when  the 
trajectory  has  gone  through  the  turning  point  and  once 
again  reaches  HI  a conditional  statement  automatically 
changes  NEND  to  one.  DE  has  been  altered  so  that  when 
NEND  equals  one  the  calculation  is  cut  off  exactly  as  if 
TOOT  had  been  reached.  Notice  that  there  is  no  fixed 
value  corresponding  to  XKIN  in  LOGDEBIV.  The  turning 
point  of  the  trajectory  serves  this  purpose  naturally. 

VEFFI  is  the  initial  value  of  the  effective 
potential.  It  can  take  one  of  three  values.  If  an  SCE/EP 
(l->2)  calculation  is  being  performed  VEFFI»V11  (Vll  is 
the  lower  potential  surface).  If  an  SCE/EP  (2->l) 
calculation  is  being  performed  VEFF1-V22  (V22  is  the  upper 
potential  surface).  If  an  SCE/ABP  calculation  is  being 


pecfornBd  VEfn-(VlltV22)/J  , The  initial  menentua,  also 
needed  for  the  calculation,  is  PMr*S0RTf2  RM(  E-VEFri ) I . 
Notice  that  the  Initial  total  energy  z is  not  always  the 
Initial  kinetic  energy.  It  depends  on  which  VEPFI  is 


the  coefficients.  The  SCE/AEP  calculation  has  the 
Itens  along  with  four  additional  coefficients. 


this  case  a dummy  value}. 

PELERR  and  ABSERR  are  the  maximum  allowable  relative 
and  absolute  errors.  DE  evaluates  expressions  for  these 

than  RELERR  or  ABSERR  the  calculation  is  either  terminated 
or  an  error  condition  will  be  indicated  in  the  output. 

The  results  reported  in  this  work  have  no  such  problems. 

A value  of  l.OE-OS  is  Che  largest  acceptable  for  results 
that  converge  consistently,  Swaminathan  and  his  covorkers 
perform  calculations  of  this  type  with  an  error  of  l.OE-07 
(Swaminathan,  19671,  Since  Che  adiabatic  equations  are 
more  demanding  l.OE-OB  was  chosen  to  give  a little  extra 
margin  of  safety.  All  that  happens  in  a practical  sense 
is  that  DE  takes  a few  extra  steps  integrating  through  the 
turning  point  region. 

lELAG  gives  Che  statue  of  the  calculation.  lELAG^l 
is  the  usual  condition  at  the  beginning  of  a new 
calculation.  The  user  is  referred  to  the  documentation 
for  Che  significance  of  the  other  values  of  IFLAG 
(shampine  and  Gordon,  1976), 

Finally,  Figures  5-10,  5-11,  and  5-12  show  the  stable 
region  of  e for  a sample  trajectory  for  each  model 
potential.  The  smaller  the  value  of  e,  the  more  nearly 
the  turning  point  region  is  truly  singular,  and  the  longer 
it  takes  OG  to  integrate  through,  very  small  values  of  t 
use  up  an  extraordinary  amount  of  computer  time.  The 
, then,  is  to  find  the  largest  value  of  e in  the 
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stable  re9ion  and  use  it.  The  figures  show  that,  at 
£■1.79  Hartree,  c«D.01  is  a reasonable  choice.  In  fact 
the  effects  of  e will  vary  with  £.  Fortunately  It  turns 
out  that  this  choice  of  t is  reasonable  for  all  the 
energies  exaisined  in  this  work  for  the  NC  and  CR 
potentials.  The  aore  reasonable  choice  for  Che  CA 
potential  at  all  energies  is  an  order  of  magnitude 
smaller , e-0 , 001 . 


All  three  graphs  show  the  analytieal 
derivative.  Equation  (5-1).  as  the  solid  line. 
The  points  (*+*)  f" 


BOTTOM  GRAPH  - 


MIDDLE  GRAPH  - 


TOP  GRAPH  - 


of  the  diabetic  potential  elements  (NO. 

All  three  graphs  show  the  analytical  second 
derivative.  Equation  (5-11,  as  the  solid  line. 
The  points  (**+)  are  obtained  numerically 
using  IHSL  Subroutines  (see  Section  5-1). 
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CHAPTEB 


HUHERICAL  RESULTS 
6-1.  Introduction 

This  chapter  is  divided  into  Couc  sections.  The 
second  section  contains  figures  showing  transition 
probability  versus  energy  produced  fron  EXACT  and  sce/ep 
(DlAB  and  ADXAB.  l->2  and  2->ll  calculatlone  for  the  NC, 
CR.  and  CA  potentials.  Acconpanylng  the  graphs  are  a few 
observations.  The  third  section  develops  the  concept  of 
the  averaged  effective  potential  and  shows  how  it  is 
applied.  The  motivation  for  using  an  average  of  effective 
potentials  is  to  include  detailed  balance 
( microreversability ) in  the  calculations.  Several  figures 
accompany  the  explanation.  The  fourth  section  shows 
results  (probability  vs.  energy)  for  the  5CE/AEP,  DlAB  and 
ADlAB,  calculations  for  the  three  model  potentials.  The 
figures  in  this  section  have  two  salient  features.  First, 
the  diabatic  and  adiabatic  curves  are  practically  Che 
same,  on  the  scale  shown  they  ace  nearly 
Indistinguishable.  This  is  the  Ideal  consequence  of  the 
fact  that  both  approaches  utilise  the  same  information. 


Second,  there  is  cemarkabla  agreement  with  the  SICACT 
results.  Detailed  balance  is  a critically  important 
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fast  oscillations 
lar^a  oaxiDuD 


In  going  fron  NC  to  CR  the  coupling, 


decreases  by  a 


factor  of  6,  but  the  CR  surfaces  cross,  the  denominator 
goes  through  aero,  the  ratio  becomes  large,  there  are 
faster  oscillations  and  a larger  maximum  In  going 

from  CR  to  CA  both  the  coupling  and  surface  separation 
diminish.  Apparently,  however,  the  ratio  is  generally 
larger  for  the  CA  case  which  oscillates  more  rapidly. 

There  is  another  trend  in  going  from  NC  to  CR  to  CA. 
The  magnitude  of  the  diabatic  coupling  decreases.  The 
smaller  the  coupling  the  more  nearly  alike  the  diabatic 
and  adiabatic  equations  become.  The  resemblance  between 
the  diabatic  and  adiabatic  curves,  even  without  detailed 
balance,  increases.  The  figures  bear  this  out.  For  the 
NC  case  agreement  between  the  SCE/EF  DIAB  (1->Z),  Figure 
6-2,  and  SCE/SP  ADlAB  (1->2I,  Figure  6-3,  curves  is  quite 
poor.  The  same  curves  in  Che  CR  case.  Figures  6-5  and 
6-6,  agree  very  well  except  at  very  low  energies.  The 
adiabatic  curves  for  Che  CA  case.  Figure  6-9,  are  included 
here  more  for  Che  sake  of  completeness  than  because  of  any 
new  information  they  night  contain.  They  are  practically 
identical  to  Che  diabatic  ones.  Figure  6-6.  Numerical 
values  show  very  small  differences  at  low  energies. 
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It  nay  also  be  observed  that  in  the  high  energy  Unit 
colliaions  take  place  too  rapidly  for  transitions  to 
occur.  This  is  evidenced  by  decreasing  transition 
probability  as  energy  increases  in  the  high  energy  regime. 
This  is  precisely  what  is  seen  in  the  CR  and  CA  figures. 
The  NC  calculation  is  not  taken  high  enough  in  energy  for 
this  phenomenon  to  be  evident.  Presumably  one  would  see 
the  same  behavior  in  this  case  also. 

The  reader  will  notice  a peculiarity  in  the  (l->2) 
calculation  for  the  CA  potential,  left-hand  graphs  of 
Figures  6-$  snd  6-9.  At  about  0.42  Hartree  the  curve 
doesn't  complete  its  cycle,  it  starts  to  descend,  but 
rises  again  well  before  reaching  sero.  As  discussed  in 
the  next  section,  the  effective  potential  for  the  SCS/EP 
calculations  generally  follows  the  potential  surface  it 
starts  out  on  to  the  turning  point  and  comes  out  displaced 
from  that  surface  in  proportion  to  the  transition 
probability.  The  lower  surface  for  the  CA  potential  has  a 
veil,  tt  would  appear  that  when  the  energy  is  increased  a 
certain  amount  the  effective  potential  has  a problem 
following  the  well.  The  fact  that  the  (2->l)  calculation, 
right-hand  graphs  of  Figures  6-B  and  6-9,  has  no 
peculiarity  supports  this  notion.  It  examines  the  same 
system,  but  doesn't  have  to  deal  with  the  well.  Since  all 
of  the  (1->21  trajectories  do  conserve  energy  the 
explanation  given  is  not  sufficient.  As  it  turns  out, 
this  monograph  supports  the 


SCE/AEP  approach 


SCE/E?  approach,  rot  that  reaaon  not  much  time  was  spent 
examining  the  phenomenon. 

The  effect  of  potential  wella  remains  a possibly 
fettlle  area  for  future  Investigation.  It  should  have 
something  to  do  with  resonances. 
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diceetion  and  ‘6  the  effective  potential  eonatcueted 
fcon  2->l  coefficients,  a single  trajectory  may  be 
calculated  from  the  averaged  effective  potential, 

^ ■ ~S~  • (6-1) 

This  trajectory  is  self-consistontly  coupled  to  two 
complete  sets  of  coefficients  Chat  are  equivalent  to  each 
other  except  that  one  starts  in  the  lower  state  while  the 
other  starts  in  the  upper.  If  the  order  were  reversed,  so 
that  the  coefficients  starting  in  the  upper  state  started 
in  the  lower  and  the  set  starting  in  the  lower  state 
started  in  the  upper,  the  results  of  the  calculation  would 
be  unaffected. 

Hence  we  postulate  that 

"To  be  successful,  the  eikonal  calculation  performed 
with  the  averaged  effective  potential  must  give  results 
such  that  “’ll"'’22’-" 

This  statement,  placed  in  quotes  for  emphasis,  can  be  used 
in  connection  with  two  mildly  different  treatments. 

In  the  first  treatment  one  performs  the  calculation, 
evaluates  Pjj  and  (Pjj  and  P22’’  confirms  that  they 
are  equal,  and  states,  therefore,  that  the  calculation  is 
successful.  This  treatment  is  used  with  the  diabetic 
eikonal  calculation,  where  the  calculation  as  it  has  been 
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described  Is  capable  of  giving  meaningful  probabilities  at 
all  R(t).  In  fact,  the  diabatic  equations  perform 
beautifully  using  this  approach.  Performing  the 
calculation  without  any  additional  constraints,  each  set 
of  coefficients  naturally  yields  probabilities  that  (to 
within  the  numerical  error  permitted  by  the  program)  equal 
each  other  at  all  internuclear  separations  and  all 
energies, 

„2 

^21  ■ '’12  ■ '’tr  ■ (6-2) 

where  refers  to  the  probability  that  the  state  of  the 

system  will  undergo  a transition  (tr  for  short).  Figure 
f-10  shows  versus  time. 

A second  treatment  has  been  developed  for  use  with 
the  adiabatic  equations.  The  approach  is  basically  the 
same  except  that  the  calculation  is  not  designed  eo  that 
^2l"^12  ^ll"^22^  B(t).  Valid  probabilities  are 

only  obtained  at  internuclear  separations  where  these 
equalities  hold.  Probabilities  at  these  locations  ace 
accepted,  all  others  are  rejected. 

why  does  the  first  treatment  work  with  the  diabatic 
equations  and  not  the  adiabatic  equations?  Notice  that 
the  diabatic  equations  for  the  coefficients  depend  only  on 
the  coefficients  A themselves  and  Equation  (4-lS)  or 
(4-18).  They  are  independent  of  whatever  effective 
potential  is  used.  The  adiabatic  equations  ace  far  more 
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explicit,  they  depend  not  only  on  A and  V®,  but  on  terma 
in  p and  u as  well.  This  leaves  a ehoioe.  One  could 
solve  for  Aj  using  Vj  and  (-  [2b  . 

using  Vj  and  Pj  (-  [2b  ) [see  Equations 

14-54)  and  (4-55)  ).  Notice  Chat  Uj.Uj  Without 

going  into  a gceat  deal  of  detail  (these  equations  becoae 
very  tcleky),  one  finds  that  this  kind  of  approach  leads 
back  in  the  direction  of  two  separate  solutions  for 
and  and  gives  pretty  much  the  sane  (inoorrect)  results 
obtained  with  the  non-averaged  effeotive  potentials. 

There  is  no  choice  but  to  seek  an  alternative. 

The  alternative  is  to  use  u and  P (from  the 
trajeotocy  equations,  the  integral  of  Equation  (4-53)  ) in 
all  of  the  adiabatic  differential  equations  for  the 
coefficients.  This  has  the  advantage  that  the  sets  of 
differential  equations  corresponding  to  the  two  sets  of 
coefficients  ace  identical  in  focn,  a situation  that  is 
known  to  work  in  the  diabetic  case,  unfortunately,  both 
sets  of  equations  are  now  incorrect.  They  each  contain  w. 
which  is  an  explicit  function  of  both  sets  of 
coefficients.  Fortunately,  this  objection  is  ceaoved  when 
the  coefficients  in  question  are  equal!  One  now  has  valid 
coefficients  that  were  calculated  using  a potential  that 
treats  collisions  starting  in  the  lower  and  upper  states 
equivalently.  For  example,  at  certain  locations. 


($-31 


In  fact,  Is  a secies  of  discrete  data  points 

cotrespondlng  to  pcobabilities  at  the  locations  where  the 
coefficients  are  equal.  It  is  a little  easier  to 
visualise  behavior  if  it  is  cast  in  the  focB  of  a 

continuous  function  whose  acceptable  values  coincide  with 
the  discrete  data  points. 

figure  $-11  shows  three  graphs  with  a conson 
horizontal  axis.  The  botton  graph  shows  both  and 
The  Biddle  graph  shows  and  Yjj.  The  top  graph  shows 
the  continuous  forn  of  The  vertical  lines  indicate 

the  locations  where  the  coefficients  and  probability  ace 
valid.  These  lie  where  the  curves  for  the  eoeftielents 
intersect.  At  sufficiently  long  tlsias  (when  the  outbound 
itajectory  is  clearly  past  the  avoided  crossing)  the 
consecutive  discreet  data  points  of  P^^  level  off  to  sone 
constant  value.  This  is  taken  as  the  post-collision 
transition  probability  at  the  specified  total  energy.  If 
one  uses  the  same  treatnent  to  separately  calculate  the 
probability  that  no  transition  will  cake  place, 
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'‘ll  * *11  *22  * *22  * *nl  * *nt  ’ ^11  * ^22  ’ ^nt  ' 

(6-4) 

and  adds  it  to  the  transition  probability,  one  does  find 


as  it  should,  while  this  approach  was  designed  to  nake 
the  adiabatic  equations  wor(c,  there  is  nothing  to  prohibit 
using  it  with  the  diabatic  equations.  Of  course  this  is 
not  preferred,  but  it  is  interesting  to  compare.  Figure 
6-12  shows  the  same  information  as  Figure  6-11,  but 
obtained  using  the  diabetic  equations.  The  two  figures 
are  very  similar,  in  this  particular  example  the  location 
of  the  post-collision  data  points  practically  coincide.  A 
notable  feature  is  that  coefficients  from  the  second  set, 
Xj2  and  V^2'  level  off  rather  than  continuing  to 
oscillate,  while  this  was  certainly  not  intended  or 
expected,  there  is  nothing  to  counter-indicate  this 
behavior . 

The  most  obvious  limitation  on  the  treatment  just 
described  is  its  inability  to  give  probabilities  that  are 
continuous  functions  of  time  (or  R(t)  ).  One  nay  attempt 
to  overcome  this  deficiency  by  connecting  Che  pea)(s  in  the 
continuous  form  of,  say,  to  get  at  least  a partial 
picture.  Figure  6-13  compares  the  continuous  probability 
from  Figure  6-10  to  peak-to-peak  curves  made  from  Figures 
6-11  and  6-12.  Starting  a little  past  the  turning  point. 
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they  ace  all  equal.  In  the  turning  point  region  itself, 
the  adiabatic  calculation  is  totally  uncellable  due  to  the 
use  of  E.  In  fact,  though  Figure  6-13  doesn't  show  it, 
the  adiabatic  peak-to-peab  curve  nay  show  erratic  behavior 
in  this  region.  The  diabatie  peak-to-peak  curve  naturally 
coincides  with  the  continuous  diabatie  probability  of 
Figure  6-10.  Since  there  is  nothing  to  be  gained  (and  a 
good  deal  to  be  lost)  perforning  diabatie  calculations 
with  Che  second  treatment,  any  diabatie  calculation 
referred  to  beyond  this  point  has  been  done  with  the  first 
approach.  If  all  that  is  required  is  the  post-collision 
probability  the  adiabatic  equations  may  be  used.  If  the 


time  evolution  of  the  probability  is  needed  t 


1 eikonal 


method  c 
representation. 


C present  only  b 


Finally,  the  actual  adiabatic  calculations  were 
performed  using  a modification  of  the  averaged  effective 
potential  in  Equation  (6-1).  Rather  than  normalising  Che 
coefficients  so  chat 
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turning  point 


the  original 
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averaged  effective  potential,  (w,  *v,  >/2,  E<juation  (6-13, 
only  function  this  «ay  with  relatively  large  e,  5.0  >e 
>0.1  . T3ieae  calculationa  behave  reasonably  in  every 
respect  except  for  the  results,  tn  this  case  the  validity 
of  the  coefficients  is  nullified  by  large  t.  on  the  other 
hand,  any  small  deviation  that  appears  in  either  of  the 
collectively  normalized  sets  at  the  turning  point  region, 
appears  to  be  simultaneously  compensated  by  an  equivalent 
deviation  in  the  other  set,  and  eventually  they  self- 
correct.  This  allows  much  smaller  values  of  e to  be  used 
and  the  coefficients  remain  valid. 

All  of  the  SCS/ACP  calculations  described  in  this 
woric  were  first  carried  out  using  the  averaged  effective 
potentials  described  above.  Fortunately  this  is  not  the 
whole  story.  The  very  fact  that  mletoreversibility  can  be 
successfully  imposed  on  the  calculations  can  be  used  to 
simplify  things  even  more.  The  averaged  effective 
potential  In  the  diabatic  representation  is 


where  w^  is  a weight  factor  assigned  to  each  level  j.  In 
the  simple  two-state  case  where  each  curve  is  given  equal 
weight  (as  above)  w^  - Wj  . 1/2  . In  general  there 
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oight  be  reasons,  such  as  level  degeneracies,  to  assign 
the  Wj  different  values. 


■ s 


■■  s- 


this  leads 


vhich  is  arrived 


♦ ^22  , (6-12) 
at  fron  the  orthonormal ity  relationship 
‘kl  - (6-13) 


which  itself  is  a consequence  of  the  foris  of  the 
differential  equations  for  the  electronic  anplltudes. 
In  the  case  being  considered 


This  little  bit  of  lueh  allows  the  trajectories  to  be 
calculated  independently  of  the  coefficients.  In  fact. 
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the  same  development  applies  to  the  adiabatic  equations  so 
that,  as  long  as  the  coefficients  ace  pcopecly  nocmaliaed, 
the  SCE/AEP  calculations  can  be  pecfocned  without  using  s! 

It  is  important  to  realize,  however,  that  when  the  w 's 
7 

are  not  equal,  the  averaged  effective  potential  will 
depend  on  the  coefficients  and  the  more  complicated 
equations  must  be  used. 


Before  moving  on  it  should  be  pointed  out  that  the 
Simplified  SCE/AEP  calculations  lake  just  as  long  as  tho. 
that  use  the  mote  general  equations.  Calculation  of  the 
coefficients,  not  the  trajectories,  is  rate  determining. 
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6-<,  SCE/AEP  Tcangltlon  Probabilities  and  Sumnary 

Figuce  £-14  shows  ths  SCE/ABP  DIAB  and  AOIAB 
transition  probabilities  versus  energy  for  the  NC 
potential.  Figures  6-15  and  6-16  show  them  for  the  CR  and 
CA  potentials  respectively.  The  diabatic  and  adiabatic 
curves  are  nearly  identical.  The  agreement  between  these 
curves  and  the  EXACT  results,  shown  here  as  dashed  lines, 


remarkably  good. 
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The  first  two  trends  support  the 
averaged  effective  potentlali  the  third  trend  supports  the 
credibility  of  the  adiabatic  approach  developed  in  this 
docunent.  By  now  there  can  be  little  question  of  either 


EXACT  calculations  ate  possible,  the  sce/ABP  calculations, 
which  appear  to  be  quite  good,  are  soaewhat  quieter.  This 
is  clearly  seen  in  the  right-hand  coluoin  in  the  above 
tables.  Pot  the  NC  and  CR  potentials  there's  about  an 
eight-fold  savings  in  tiae  with  the  DlAB  ealeulation  and 

calculation,  in  the  worst  ease  in  this  regard,  the  CA 
potential  , there's  still  a three-and-one-half-fold 
savings  with  DIAB  and  two-fold  with  ABIAB.  It  is  worth 
knowing  that  these  programs  were  written  solely  to  develop 
and  perfect  the  calculations.  Even  a rudimentary  attempt 
at  optimising  the  code  would  result  in  some  greater 
savings. 

The  tables  reveal  two  more  interesting  facets  of  the 
AEP  calculations.  First  of  all,  the  SCB/AEP  ADIAB 
calculations  invariably  take  longer  than  the  SCE/AEP  DIAB 
calculaticns.  This  adds  weight  to  the  previously 
discussed  notion  that  the  adiabatic  approach  should  only 


• The  CA  calculations  is  slower  than  the  HC  and  CB  because 
the  reduced  mass  is  30  tines  heavier.  Also  the  c used  in 
the  A01A8  CA  calculation  is  one-tenth  the  value  used  in 


be  applied  when  Che  dlabaCic  approach  is  unseCisfactocy. 
John  Tully  suggesCs  chat  there  ace  cecCain  low  energy 
situations  where  this  night  be  the  case  (Tully,  1980). 

The  second  facet  is  a little  more  curious,  Che  AEP 
calculations  take  very  little  more  (and  in  many  cases 
actually  lese)  computer  time  than  the  EP  calculations  even 
though  nearly  twice  as  many  differential  equations  ace 
being  solved.  The  design  of  DB  is  such  that  the  calls  to 
and  from  the  various  subroutines  take  proportionately 
longer  compared  to  the  time  spent  within  those  routines. 
This  means  doubling  the  number  of  differential  equations 
will  not  generally  double  the  computet  time.  But  this 
doesn't  fully  explain  the  diminutive  nature  of  the 
increases  seen,  and  it  doesn’t  explain  the  decreases  at 
alll  These  appear  to  be  brought  about  by  the  increased 
symmetry  of  the  SCE/ABP  equations.  In  fact,  in  cases 
where  there  is  a decrease  in  CPU  time,  energy  conservation 
is  actually  several  orders  of  magnitude  better  than 
dictated  by  error  tolerance!  Energy  conservation  is 
obtained  from  Che  simple  expression 


where  E is  the  initial  total  energy,  pVim  is  the  current 
kinetic  energy,  and  u is  the  current  potential  energy.  In 
the  usual  case  when  the  error  parameters  input  to  DB  are 
set  to  10  , a will  eventually  settle  to  10“*  or  possibly 

10  In  cases  where  the  increased  number  of  equations 
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pioduees  a reduction  In  conputer  tioe  a will  settle  to 
anywhere  from  ID  to  10  without  any  added  incentive 
in  the  ptogran.  D£  uses  polynomial  apptoalmatlons  to 
predict  the  behavior  of  the  equations  it  then  solves 
numerically.  Apparently  there  is  some  matching  of  the 
polynomial  predictions  when  two  sets  of  coefficients  ace 
used  that  makes  it  easier  for  DE  to  handle  the  whole  set 
of  equations. 

Overall,  Che  SCE/AEP  calculations  are  more  accurate 
than  the  SCE/EP  calculations  by  a vast  margin.  They  are 
also,  as  indicated  in  the  last  paragraph,  cost  effective 
in  terns  of  computer  time,  producing  a slight  increase  in 
some  cases  and  an  actual  savings  In  others.  The  present 
work  definitely  encourages  the  use  of  averaged  effective 


potentials. 


CONCLUSION 


After  a 9eneral  introduction  to  the  eikonal  method 
(Hicha,  1983),  the  formaliein  for  previously  performed 
three-dinenaional  elaetic  scattering  calculations  was 
presented  along  with  a few  sample  results  (Cohen,  1988). 
Following  that  was  the  development  of  the  ei)conal 
formalism  for  inelastic  scattering  in  one  variable.  A )rey 
part  of  Che  effort  was  aimed  at  applying  the  method  in 
both  Che  diabetic  and  adiabatic  representations.  The 
presentation  made  use  of  three  model  two-state  potentials 
selected  Co  cover  a range  of  conditions.  There  was  a non- 
crossing potential  (NO,  a potential  with  curve  crossing 
in  the  repulsive  region  (CH),  and  a potential  with 
crossing  in  the  attractive  region  (CA).  Development  of 
Che  diabaCic  equations  followed  that  of  Hicha  and 
Swaminathan  (Hicha  and  Swaminathan,  1986).  The  adiabatic 
equations  presented  in  this  wor)c  were  formed  by  applying  a 
unitary  transformation  to  the  diabaCic  equations.  They 
displayed  some  very  interesting  features,  both  formally 
and  in  their  computational  characteristics. 

The  most  obvious  difference  in  Che  adiabatic 
equations  is  Che  addition  of  several  terms  to 
schroedinger ' s equation  for  nuclear  motion  that  result 
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£com  the  unitaty  transfocnation.  When  the  dlahatic 
equation,  already  simpler  than  the  adiabatic,  ie  switched 
to  the  time  dependent  form,  it  can  be  simplified  further 
yet  by  an  appropriate  chanqe  in  the  phase  of  the  eikonal 
wavefunctlon.  This  amounts  to  a rearrangement  so  that 
unwanted  espressions  become  part  of  the  exponential 
factor,  which  is  then  divided  out.  While  the  sane 
strategy  ie  used  to  simplify  the  adiabatic  equation 
somewhat,  it  cannot  be  used  to  nearly  the  same  extant.  A 
leas  obvious  difference  stems  from  the  appearance  of  a 
momentum  variable.  P,  in  the  denominator  of  the  adiabatic 
equations  for  the  coefficients.  This  leads  to  two 
separata  singularities  in  Che  vicinity  of  the  turning 
point,  where  P goes  to  aero.  These  are  avoided  by  the 
addition  of  a small  imaginary  increment,  e,  to  P.  if  c is 
too  small  the  calculation  will  take  too  long,  if  e is  too 
large  the  eoeffieients  will  not  be  valid.  Testa  wore  made 
for  the  appropriate  choice  of  s. 

Other  differences  concern  calculation  of  the 
trajectories.  Hamilton's  equation  for  the  temporal 
derivative  of  the  momentum  requires  the  spatial  derivative 
of  (in  this  easel  the  effective  potential.  It  might  be 
expected  that,  since  the  effective  potential  contains  the 
coefficients,  this  would  require  the  spatial  derivatives 
of  the  coefficients.  In  the  diabetic  case,  however,  the 
coefficients  can  also  be  arrived  at  using  Hamilton's 
equations.  Because  of  the  symmetric  relationships  among 
coefficients  formed  in  this  way  (see  Equation  (4-18)  ), 
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the  spatial  daclvatives  that  appear  in  3u*^/9B  cancel  one 
another  (see  Equation  (4-2E)  ).  in  the  adiabatic  case,  on 
the  other  hand,  the  spatial  derivatives  o£  the 
coefficients  are  needed.  Also,  the  adiabatic  effective 
potential,  even  after  application  of  the  shoct-vavelength 
approximation,  contains  a term  in  P.  This  momentum 
dependent  potential  requires  a little  special  attention. 

Finally,  the  time-dependent  adiabatic  equations 
contain  certain  non-Hermltian  terns.  This  means  that  the 
expression  used  for  probability,  a'^a.  Is  not  automatically 
conserved.  Normaliration  must  be  provided  to  ensure  that 


All  of  these  oddities  and  complications  certainly 
make  the  adiabatic  formalism  more  difficult  as  well  as 

more  demandinq  equations  are  less  hardy  and  resilient  than 
the  simple  symmetric  diabetic  equations,  and  most 
importantly,  they  take  more  computer  timel  Since  the  same 
information  goes  into  both  the  SCE  DIAB  and  APIAB 
calculations.  Ideally,  the  same  results  ace  expected  from 
each.  In  fact  (as  something  of  a triumph)  this  is  true 
for  the  calculations  presented  in  this  work, 
unfortunately  however,  the  adiabatic  calculations  take 
longer  than  the  diabatic  calculations  for  every  single 
example.  All  of  the  reasons  In  this  paragraph  combine  to 
make  the  adiabatic  approach  generally  less  desirable  than 
the  diabatic  one.  One  suspects,  of  course,  that,  with  the 
infinite  variety  of  chemical  systems  and 


potentials 


dealt  with,  sone  will  be  note  anenable  to  the  adiabatic 
tceatstant. 

In  any  caae,  the  successful  development  and 
application  of  the  adiabatic  eikonal  equations  presented 
in  this  work  adds  to  the  dynamielst' s understanding  of  his 
subject.  This  author,  at  least,  envisions  the  possibility 
of  applying  the  eikonal  method  to  equations  not 
transformed  from  diabatic  ones.  Since,  as  John  Tully 
points  out  (John  Tully,  1976),  curve  crossings  are  only 
special  one-dimensional  cases  of  seams  where  adiabatic 
surfaces  approach,  this  might  prove  very  interesting  in 
multidimensional  problems.  Alternatively,  one  can  always 
apply  the  methods  of  Olson  and  nicha  (Olson  and  nicha, 
1982)  to  develop  an  appropriate  diabatic  potential. 

once  the  usual  eikonal  calculations  were  successfully 
carried  out  in  both  representations,  the  question  of 
improvements  came  up.  Since  the  excitation  probabilities 
did  not  match  Che  de-excitation  probabilities,  as  they  do 
in  exact  calculations,  an  averaged  effective  potential  was 
introduced  to  impose  the  condition  of  raicrocever sibllity" , 
The  effective  potentials  for  Che  excitation  and  de- 
excitation processes  were  calculated  simultaneously,  along 
with  both  sets  of  coefficients,  and  trajectories  were 
constructed  from  their  average.  The  improved  agreement 
with  exact  calculations  was  remarkable  (see  figures  6-14, 

• A more  rigorous  treatment  of  microreversibility,  using 
variational  methods,  was  developed  by  nicha  and  Gasdv 
IHicha  and  Gasdy,  1967). 
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6-15,  and  6-16}.  Calculations  that  war«  only 
qualitatively  correct  using  an  effective  potential  becane 
seDi-quantitatively  to  quantitatively  correct  using  the 
averaged  effective  potential.  The  added  cost  in  conputer 
tine  of  using  Che  averaged  effective  potential  was  ninioal 
to  nonexistsnt  due  Co  an  increase  in  the  symnecry  of  the 
equations. 

A mildly  surprising  discovery  was  that,  due  to  the 
orthonormal ity  of  the  electronic  functions,  the  averaged 
effective  potential  for  the  case  where  the  curves  are 
given  equal  weight  is  simply  the  average  of  the  surfaces. 
In  this  particular  case  the  coefficients  don’t  appear  in 
the  averaged  effective  poCentlall  since  these 
trajectories  can  be  calculated  independently  of  the 
coefficients,  Che  adiabatic  calculations  for  this  case  no 
longer  require  c!  sadly,  calculations  performed  this  way 

coefficients  in  the  averaged  effective  potential.  The 
tine  is  taken  up  in  calculating  the  coefficients 
themselves.  In  Che  general  case  reasons  might  be  found  to 
assign  the  curves  different  weights,  leaving  no  choice  but 
to  use  the  coefficient  dependent  potential. 

The  improved  results  at  little  or  no  extra  cost,  along 
with  the  information  from  the  preceding  paragraphs,  lead 
Co  Che  conclusion  that  Che  method  of  choice  uses  the 
averaged  effective  potential  In  the  diabatlc 
representation  (SCE/AEP  OIAB). 


The  strength  of  the  eikonel  method  lies  not  in  one- 
variable  calculations,  which  can  be  performed  exactly,  but 
in  oultl-vatlable  calculations,  which  cannot.  Clifford 
Stodden  has  used  the  method  to  treat  the  photodissociation 
of  methyl  iodide  (an  interesting  application  in  itself)  to 
two  variables  (Stodden.  19S7), 

Extending  the  method  is  the  heart  of  simplicity.  A 
set  of  coefficients  is  added  for  each  new  surface.  The 
total  probability,  A*A.  is  still  set  to  unity  and 
conserved.  The  sum  of  the  squares  for  each  set  of 
coefficients.  aJAj  - +y?,  represents  the  probability 

that  the  system  is  on  that  surface,  i.  An  averaged 
effective  potential  might  easily  be  formed  from  three  or 
more  surfaces. 

t least  conceptually,  extension  of 
e inelastic  scattering  information  in 
lis  would  give  the  scattering  angle  as 
To  accomplish  this  one  would  use  the 
the  Internuclear  separation  R and  the 
the  momentum  variable  P.  Since  the 
»enta  would  not  be  equal,  there  would, 
ise,  be  a range  of  final  momenta  for 
Also  the  trajectories  would  not 
remain  in  a single  plane.  Although  these  complications 
would  add  to  the  complexity  and  cost,  in  computer  time,  of 
the  calculations,  the  procedure  remains  eminently 
practical.  A largo  part  of  the  challenge  would  lie  in 
constructing  appropriate  potential  surfaces.  Central 


Another  simple, 
the  method  would  gi' 
three  dimensions. 

initial  and  final  me 
unlilte  the  elastic  c 
each  scattering  angle. 
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poCentials  obtained  by  averaging  over  scattering 
experiments  might  be  used  for  model  calculations. 

In  19B4  Olson  and  Garrison  performed  preliminary 
elkonal  calculations  for  the  scattering  of  Na  from  W(IIO) 
surfaces  (Olson  and  Garrison,  1985).  In  1987  more 
sophisticated  model  calculations,  taking  surface  motion 
into  account,  were  applied  to  collisions  of  gas  atoms  with 
insulator  surfaces  < Swaminathan , Garrett,  and  Rurthy, 
1988).  Equations  of  motion  characterising  the  surface 
were  formed  from  a stochastic  ceduction  of  the  generalised 
Langevln  equation  (GLE).  for  systems  studied  there  was  an 
appreciable  loss  of  gas  atom  energy  to  the  surface, 
clearly  related  to  the  inelastic  transition  probability. 
Possible  extensions  of  Che  theory  Include  the  treatment  of 
photochemical  processes  on  surfaces,  and  reactions  with 
and  on  surfaces.  In  1988  roCationally  electronically 
inelastic  scattering  of  HO  from  Ag(lll)  was  treated  using 
Che  eikonal  method  and  results  compared  very  well  to  exact 
calculations  (Rice,  Garrett,  and  Swamlnathan,  1989). 

These  calculations  seem  to  be  limited  only  by  the 
availability  of  potentials  along  with,  as  with  the  GLE, 

Che  various  other  parameters  Chat  are  used.  For  the 
moment  at  least,  that  is  where  things  stand.  While  our 
handling  of  it  can  always  be  improved,  the  basic  eilconal 
approach  is  wholesome  and  theoretically  sound  and,  as  long 
as  correct  input  is  available,  useful  results  will  he 


generated. 


APPENDIX  A 
KEY  TO  PROGRAM  NAMES 


APPENDIX  B.  EIXONAL 


PROCBAn 


C THIS  PROGRAK  CALCULATES  ELECTRONICALLY  INELASTIC 
C TRANSITION  PROBABILITIES  FOB  TBE  RC  POTENTIAL  (CROSSING 
C IN  REPULSIVE  REGION)  USING  THE  SELF  CONSISTENT  EIXONAL 
C APPROACB  IN  THE  DIABATIC  REPRESENTATION  WITH  AN  AVERAGED 
C EFFECTIVE  POTENTIAL.  IT  IS  PART  OF  "COMPUTATIONAL 
C STUDIES  OF  ELECTRONICALLY  DIABATIC  MOLECULAR  COLLISIONS" 

C A DOCTORAL  DISSERTATION  SUBMITTED  BY  JOEL  COHEN  TO  TBE  ' 
C UNIVERSITY  OP  FLORIDA  IN  OCTOBER  1990. 

C SEE  CHAPTER  II  FOR  THE  GENERAL  EIXONAL  EQUATIONS,  SECTION 
C 4-2  FOR  THE  DIABATIC  EIXONAL  EQUATIONS,  TABLE  4-1  AND 
C FIGURE  4-1  FOR  THE  NC  POTENTIAL,  AND  SECTION  5-3  FOR  AN 
C EXPLANATION  OF  THE  AVERAGED  EFFECTIVE  POTENTIAL 
C THE  AVERAGED  EFFECTIVE  POTENTIAL  IS  ONE  HALF 
C TBE  SUM  OF  THE  EXCITATION  AND  DE-EXCITATION  EFFECTIVE 
C POTENTIALS.  FUNCTIONS  TO  BE  EVALUATED  AT  POINTS  ALONG  THE 
C TMJECTORY  ARE  STORED  IN  Y(I).  THE  INITIAL  VALUES  ARE 
C FOUND  IN  THE  MAIN  PROGRAM.  A DIFFERENTIAL  EQUATION  SOLVER 
C (DE)  IS  CALLED  TO  SIMULTANEOUSLY  GENERATE  TRAJECTORIES 
C AND  TRANSITION  AMPLITUDES  STARTING  AT  THE  INITIAL 
~ ‘■ISTED  IN  MAIN.  OE  USES  SUBROUTINE  FPRIMB 

C Of  yf j^'^”^****®  TP(I),  THE  ANALYTICAL  TEMPORAL  DERIVATIVES 
C YU)  IS  THE  INTERNUCLEAR  SEPARATION. 

0 Y(2)  IS  THE  RELATIVE  MOMENTUM. 

C YO)-Y(S)  ARE  THE  COEFFICIENTS  OF  EXCITATION. 

C Y(7)-Y(101  ARE  THE  COEFFICIENTS  OF  DE-EXCITATION. 

C IT  IS  NOT  NECESSARY  TO  REPRODUCE  THE  CODE  FOR  OE, 

C THAT  IS  AVAILABLE  PROM  OTHER  SOURCES.  ONLY  ORIGINAL 
C MATERIAL  REQUIRED  BY  THE  PRESENT  WORK  IS  LISTED  BELOW. 


:C  MAIN  PROGRAM 
: INSERT  INITIAL  CON 
: TO  BE  EVALUATED  BY 
IMPLICIT  REAL'S 


EXTERNAL  F 
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r DERIVATIVE  A 


DHllDR— ALPHAll'Hll 
H221-A23«DEXF<-ALPHA22o<  R) ] 

H22-H221-fH222W20 

DB22DB— ALPRA22*H221-BETA22*H222 
TANA-TAHH(ALPHA12*(R-B12 ) ) 


C ATHA/AtA  IS  THE  EFFECTIVE  POTENTIAL  FOB  EXCITATION 


C ATHAS/ATAS  IS  THE  EFFECTIVE  POTENTIAL  FOB  DE-EXCITATION 


C VBFF  IS  THE  AVERAGED  EFFECTIVE  POTENTIAL 
VEFF-(ATHA/ATA)*(ATHAS/ATAS)  )/2. 
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